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FOREWORD 


This report was prepared by the McDonnell Aircraft Coapany (MCAIR), a 
division of the McDonnell Douglas Corporation, St. Louis, Missouri for the 
Rational Aeronautics and Space Adninistration, Dryden Flight Research Center, 
Edwards, California. The study was performed under RASA Contract NAS4-2364, 
"F— 15 Inlet/Engine Test Techniques and Distortion Methodologies Study." 

The work was performed from March 1977 through February 1978 with Mr. Jack 
Nugent (NASA/ Dryden) as Program Monitor and Mr. Harvey Neumann (NASA/Levis) 
as Technical Monitor. Special acknowledgement is d e Mr. T. Putnam (NASA/ 
Dryden) for his constructive criticisms and suggestions. 

The effort at McDonnell Aircraft Coapany was conducted under the tech- 
nical leadership of the Engineering Technology Division. In addition to the 
authors listed on the cover, other MCAIR personnel that made significant 
contributions to this program were Mr. Edward Smith, Mr. Lee We laser and 
Mr. Mark Sawyer. Special acknowledgement is due Mr. Hershel Sans for his 
reviews and suggestions. 

Significant subcontract support was provided by Mr. Uayne Walter and 
Mr. Lew Hayward of Pratt & Whitney Aircraft (P&WA) , Government Products 
Division, under the direction of Mr. Frank Thompson. 

This report consists of nine volumes. Technical discussions or the 
program, results and Appendices A and B are presented in Volme I (NASA CR 
144866) . Appendices C through J are presented in Volume II through IX (NASA 
CR 144867-144874) which present the distortion analysis plots and the assoc- 
iated statistical functions used for the analyses. 
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SUMMARY 


Recent emphasis on increased maneuverability r. .rements for fighter 
aircraft has necessitated an extensive engineering development effort be 
directed towards inlet/engine compatibility. Inlet/engine compatibility must 
be assessed early in the aircraft development program to allow necessary 
inlet and engine design modifications to be defined and implemented at minimum 
cost impact. This early assessment of ir.let/cngine compatibility is determin- 
ed by engine stability audits computed using inlet distortion levels trom 
subscale inlet model data and engine sensitivities to inlet distortion. 
Therefore, the accuracy with which subscale inlet model distortion levels 
predict flight test vehicle distortion levels is a crucial element in assess- 
ing inlet /engine compatibility. 

The primary goal of this distortion methodologies study was to determine 
if time variant distortion data taken from a subscale inlet model can pre- 
dict peak distortion levels for a full scale flight test vehicle. The data 
base used to accomplish this goal was collected in separate programs by MCAIR 
and NASA/Dryden. Subscale and full scale wind tunnel data were collected by 
MCAIR during the F-15 development program, and flight test data were collected 
by NASA/Dryden during the >.'ASA F-15 inlet/engine compatibility flight test 
program. This data base has a Mad, number range of 0.4 to 2.5 and an angle 
of attack range from -10 degrees to +12 degrees. 

The primary objectives accomplished in meeting the overall program goal 
wer_ to determine the effects on peak distortion of: (1) Reynolds Number/ 

scale, (2) engine presence and (3) frequency content. In addition, the ''apa- 
bility of the P&WA stability audit system to predict engine stalls was 
evaluated, and the capability of Melick’s procedure. Reference (1), to pre- 
dict peak time variant distortion levels was evaluated. Using the Pratt and 
Whitney Aircraft distortion descriptor, Ka 2 > the data indicate the following 
significant results for the F-15/F100 inlet/engine propulsion system. 

o Peak time variant distortion from |^bscale inlet model wind tunnel 
tests are representative of full scale flight test distortion. 

o The time variant pressure data of this study are random stationary 
data, thereby allowing valid statistical analyses to be conducted. 

o The effect of the engine presence on total pressure recovery, peak 
time variant distortion and turbulence level is small but favorable. 

o The Reynolds number/ scale evaluation indicates a general trend of 
increasing total pressure recovery, decreasing peak time variant fan 
distortion and decreasing turbulence with increasing Reynolds number/ 
scale. 

o The frequency content evaluation indicates that peak time variant 
fan distortion and turbulence increase with increasing filter cutoff 
frequency for all of the data evaluated in this study. 

o The capability of the Pratt & Whitney Aircraft stability audit system 
to Dredict engine stalls has been verified for both stall and non- 
stall flight test conditions. 
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c Predictions of peak distortion values using Melick's procedure are 
accurate to 11.3 percent average error for fourteen data points 
having nominal turbulence levels and are accurate to 20 percent 
average error (the maxiiMn error approaches 40 percent) for eight 
data points having high turbulence levels. 
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HfTRODUCT IOH 


Recent eaphasis on increased maneuver ability requirements for fighter 
aircraft has necessitated an extensive engineering development effort be 
directed towards inlet/engine compatibility. Inlet/engine compatibility 
must be assessed early in the aircraft development program to allow neces- 
sary inlet and engine design modifications to be definea and implemented at 
minimum cost impact. This early assessment of inlet/engine compatibility is 
determined by engine stability audits computed using inlet distortion levels 
from subscale inlet model data and engine sensitivities to inlet distortion. 
....'•fore, the accuracy with which subscale inlet model distortion levels 
predict flight test vehicle distortion levels is a crucial element in asses- 
sing inlet/ engine compatibility. 

During the development program of the F-15 air superiority fighter air- 
craft, the time variant pressure data were analyzed as required to ensure 
the F-15/F10C to be a compatible inlet/engine system prior to fleet intro- 
duction, but no in-depth analyses were conducted to investigate distortion 
methodologies and test techniques. This distortion methodologies and test 
techniques study utilizes the extensive MCAIR F-15 development test data base 
and the data base collected during the HASA/ Dry den F-15 inlet, L.. 0 »ae com- 
patibility flight test program. The wind tunnel data base was collected 
during the period of March 1970 through March 1972 using a subscale inlet 
model (l/6th Scale) , a full scale inlet model without engine (FSCP) and a 
full scale inlet model with engine (FSE) which were all tested at the Arnold 
Engineering Development Center (AEDC) at Tullahoma, Tennessee. The wind 
tunnel test conditions simulated both level flight and high energy maneuver- 
ability flight conditions. The flight test data base was collected in the 
NASA inlet/ engine compatibility flight test program during the period of 
April 1976 through June 1976 using F-15 aircraft number 71-0281 which was 
tested at the KASA/Drvden Flight Research Center at Edwards Air Force Base, 
California. The aircraft test conditions duplicated the selected wind tunnel 
test conditions for both high energy maneuverability and level fli^it test 
conditions. These two data bases were integrated to provide the data matrix 
used in this study. 

The F-15 inlet system is a two-dimensional, external compression, over- 
head three rang* variable capture area design with an internal diffuser ramp. 
This unique design is capable of rotating the entire upper cowl and compres- 
sion ramp system of the inlet , thus providing a variable capture area and 
aligning the inlet with the approaching airflow over a wide range of Mach 
numbers, angles of attack, and engine airflows. The Pratt and Whitney Air- 
craft (P&WA) FI 00 turbofan engine, is a twin spool, high pressure ratio, 
lightweight engine with mixed flow augmentation and variable area convergent- 
divergent exhaust nozzle, the fan moduli has three stages and the compressor 
has ten stages. The fan and compressor are each driven by dual turbine 
stages and the fan and core airflows are mixed in the augmentor prior to 
afterburning. 

The primary goal of this study was to determine if subscale inlet model 
peak time variant distortion data will accurately predict full scale flight 
test peak time variant distortion data. To achieve this goal, several speci- 
fic technical objectives were established including: (1) evaluate the 
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quality of the tiae variant data, (2) determine the effect of engine pre- 
sence on distortion levels, (3) establish Reynolds number/scale effects on 
distortion levels, and (4) investigate the effect of frequency content on 
distortion levels. 

The first sections of this report provide descriptions of the inlet 
models, test programs and the associated data acquisition and data reduction 
systems. The following sections describe the evaluation of the time variant 
data quality, effect of (1) engine presence, (2) Reynolds ntasber/scale and 
(3) frequency content. In addition, the capability of the P&WA stability 
audit system to predict engine stalls is evaluated and the capability of 
Melick's procedure. Reference (1), to predict most probable peak time variant 
distortion values are evaluated. The appendices contain a sample calculation 
of the distortion descriptor, K* 2 » the distortion analysis plots and the 
associated statistical functions used for the analyses. 
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MODEL DESCRIPTIONS 


The tine variant pressure data from which the peak values of time var- 
iant distortion were calculated for this study were collected using a l/6th 
scale inlet wind tunnel model, a full scale inlet wind tunnel model tested 
with and without an engine and the full scale flight test vehicle, aircraft 
71-0281. Schematics of these three test articles are shown in Figure 1. The 
applicable test series, test dates and test facilities arc as follows: 


1 

Inlet 

Model 

HCAIR 

Test 

Series 

Test 

Date 

Test 

Facility 

Facility 
Test Series 

Aircraft No. 
71-0281 

FLT 

June 1976 

NASA/DRYDEN 
Edwards, Calif. 

NASA F-15 Inlet/ 
Engine Compatibility 
Test 

Full Scale 
With Engine 

FSE 

Mar 1972 

AEDC* 

Tullahoma, Tenn. 

TF269/SF144 

Full Scale 
Cold Pipe 
(Without Engine) 

FSCP 

Nov 1971 

AEDC* 

Tullahoma, Tenn. 

SF144 

l/6th Scale 

VIII 

June 1971 

AEDC* 

Tullahoma, Tenn. 

SF139 

l/6th Scale 

VII 

July 1970 

AEDC* 

Tullahoma, Tenn. 

TF242/SF130 


*Arnold Engineering Development Center 


The general physical characteristics of the F-15 inlet configuration, 
full scale flight test vehicle, full scale inlet wind tunnel model, l/6th 
scale inlet wind tunnel model, and engine face total pressure rake are dis- 
cussed below. 


F-15 Inlet Configuration 

The inlet configuration is a two-dimensional, external compression, 
overhead three ramp variable capture area design with an internal diffuser 
ramp, and the raap and bypass system geometries are illustrated in Figure 2. 
The inlet consists of the three forward compression ramps and a diffuser or 
fourth ramp. The cowl rotation system rotates the first compression ramp 
assembly about a pivot point near the cowl lip. The second, third and 
fourth ramp systems are slaved together and are positioned simultaneously. 

The boundary layer bleed systems consist of the second and third ramp porous 
surfaces, throat slot bleed and inboard and outboard sidcplate bleeds. Hie 
throat slot in conjunction with the bypass door acts as the inlet bypass sys- 
tem for inlet/engine matching. 
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Full Scale Flight Test Vehicle 


The flight test data for this study was obtained using F-15 aircraft 
number AF 71-0281 on bail to NASA from the U.S. Air Force and is shown in 
Figure 3. This aircraft which is fully instrumented was used initially as 
the Flight Test Propulsion Subsystem Development Aircraft for the MCAIR 
F-15 flight test development program. 

On the F-15 aircraft, the second, third and fourth (diffuser) inlet 
ramps are linked together and are positioned by a utility hydraulics power 
actuator. The variable capture area feature is accomplished by a single 
actuator which rotates the cowl (including the first, second and third ramps) 
about the lower cowl lip. The bypass door is positioned by a third hydraulic 
actuator. 

The positioning of the three actuators is accomplished by an Air Inlet 
Controller (AIC) which uses freestream Mach number, inlet throat Mach number, 
total temperature, and angle of attack signals from the aircraft sensors. 

In addition, this aircraft is equipped with manual air inlet controller 
capability (not a production item) to allow inflight pilot control of the 
ramp and bypass door positions. 

The F-15 aircraft is powered by two Pratt and Whitney Aircraft F100-PW- 
100 turbofan engines. The F100 engine is a twin spool, high pressure ratio, 
lightweight engine with mixed flow augmentation and a variable area conver- 
gent-divergent exhaust nozzle. The fan module has three stages and the com- 
pressor has ten stages. The fan and compressor are each driven by dual tur- 
bine stages and the fan and core airflows are mixed in the augmentor prior 
to afterburning. The combustor is an annular design. The F100 engines used 
for this study were Serial Numbers P690063 and P690081. Engine S/N 063 
which is a pre-production engine and denoted as an F100(2 7/S) was retrofitted 
with a production fan and engine controls. This engine was installed in the 
left engine nacelle and instrumented with a high response engine face rake. 

Full Scale Inlet Wind Tunnel Model 

The full scale inlet wind tunnel model, shown in Figure 4, was tested 
with and without an engine installed, therefore providing the data necessary 
to evaluace the effect of the engine presence on peak time variant distortion 
levels. This full scale inlet test article includes the basic components of 
the inlet, inlet controls, engine and transonic and supersonic partial fore- 
bodies. 

The full scale inlet model, including the entire diffuser section, repre- 
sents the true aerodynamic configuration of the aircraft. This Inlet model 
is fabricated of flightweight hardware and the aircraft inlet description 
of the previous section is an accurate description of this inlet model. 

The inlet ramps, variable cowl and bypass door were positioned using the 
same model AIC as on the flight test aircraft. The AIC used the freestream 
Mach number, inlet throat Mach number, total temperature and angle of attack 
signals from either the model or wind tunnel sensors to position the ramps, 
variable cowl and bypass door. This AIC also had a manual mode (same as the 
flight test AIC) which allowed the test operator to vary the inlet ramps, 
variable cowl and bypass door while the tunnel was operating. 
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For the full scale inlet testing with engine, Pratt and Whitney Aircraft 
(P&WA) furnished the XF100-PW-100 ground test engine designated as XD-13. 

This engine is of the same basic design as the prototype and production F100 
engines with minor differences in the aerodynamic gas path. For the full 
scale inlet testing without engine , a flow control /measurement system, de- 
noted as the cold pipe, was used in place of the engine, i.e. full scale 
cold pipe. The size, mounting and face geometry duplicated the P&WA engine. 

The complete aircraft forebody could not be simulated at full scale be- 
cause of tunnel size limitations. Therefore, transonic and supersonic partial 
forebodies were designed to simulate the local flovfield at the inlet. Pre- 
vious inlet development testing using the l/6th scale model indicated that 
transonically over the range of angle of attack and zero sideslip covered 
in the full scale wind tunnel testin', the inlet performance and distortion 
characteristics were not significantly affected by the presence of the for- 
ward fuselage. Therefore, a small partial forebody was designed, which was 
primarily required to function as a fairing forward of the fuselage. The 
supersonic partial forebody shape was defined during the F-15 Forebody Flow- 
field Development Test. Flowfield comparisons between the full fore'oody and 
the partial forebody showed that the full forebody generated more downwash 
at the inlet plane. This downwash was duplicated during the full scale inlet 
testing with the partial forebody by rotating the entire model negatively, 
with the amount of rotation as a function of fuselage angle of attack. The 
full scale test article and supersonic partial forebody are shown in Figure 4. 

The boundary layer bleed system was identical to that on aircraft 71-0281. 
However, a test peculiar auxiliary bypass door was built into the test article 
forward of the engine face. The auxiliary door eliminated inlet instabilities 
during engine startup and windmilling. During normal engine operation, the 
auxiliary door closed to make a smooth diffuser line. All bleed systems on 
the full scale model were individually calibrated prior to testing. 

l/6th Scale Inlet Wind Tunnel Model 

The l/6th inlet wind tunnel model was used during the F-15 inlet develop- 
ment program. The model consisted of a fuselage forebodv assembly, a remote- 
ly controlled left inlet and an inlet to engine duct that is internally 
simulated to the engine compressor face station. The main duct aft of the 
engine face station is equipped with a remotely actuated mass flow plug. 

The tunnel installation of this model is shown in Figure 5. 

The l/6th scale inlet model, including the entire diffuser section, 
represents the true aerodynamic configuration of the aircraft. The cowl 
rotation system, which includes the first ramp assembly, and the second, 
third and fourth ramp system is identical to Che full scale ramp system and 
can be remotely positioned using two independent drive systems. The inlet 
assembly can be tested isolated or with the forebody. 

The aerodynamic configuration of the boundary layer breed system is 
identical to the full scale test articles, however, bleed flow ducts were 
attached to the bleed exits and these ducts were fitted with remotely ac- 
tuated mass flow plugs. The ducts could be removed if desired. 
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Engine Fa e Total Pressure Rake 


The engine face rake configuration has eight legs and six rings equal 
area weighted with the high response and steady state pressure probes located 
side by side. The l/6th scale configuration is an exact scale version of 
full scale engine face rakes as shown in Figure 6. The full scale engine face 
rake was used for the full scale inlet wind tunnel testing and also on the 
flight test vehicle. 
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TEST DESCRIPTIONS 


Valid comparisons of peak time variant distortion data from flight test 
and wind tunnel test require that the inlet configurations, test conditions, 
analysis times and inlet airflows, be equivalent. The inlet configurations, 
described above, are identical for each of the test articles. The techniques 
used to obtain comparable test conditions, analysis times and inlet airflows 
of flight and wind tunnel test articles are discussed below. 

Flight Test Procedure 

The data base for the l/6th scale inlet and the full scale inlet wind 
tunnel models existed prior to the NASA F-15 Inlet/Engine Compatibility Test 
Program. Test conditions were selected from these data that would provide a 
high level of turbulence and steady state distortion. Some of these selected 
test conditions have extreme aircraft attitudes such as -10 degrees angle of 
attack and +10 degrees angle of sideslip at 0.9 Mach number and -4 degrees 
angle of attack at 1.6 Mach number. While wind tunr.el models can be posi- 
t:' led in these extreme attitudes for indefinite periods of time, the aircraft 
cannot fly at these flight conditions indefinitely and still maintain Mach 
number and altitude. To achieve these extreme attitudes, the aircraft had 
to be maneuvered into the desired test condition. For the -10 degree angle 
of attack and +10 degree angle of sideslip at 0.9 Mach number. Figure 7, 
the aircraft was first inverted and then sideslipped. For the 1.6 Mach 
number at -4 degrees angle of attack, the aircraft was accelerated to the 
desired Mach nurber, then "pulled-up" to positive angle of attack (increasing 
altitude and decreasing Mach number) and then "pushed over" (decreasing alti- 
tude, increasing Mach number) to achieve -4 degrees angle of attack, 1.6 
Mach number at the correct altitude as shown in Figure 8. 

Analysis Time 

The available analysis times (time at desired flight condition) for the 
extreme flight maneuvers are relatively short since simultaneous matching of 
the first ramp, variable cowl position, second, third and fourth ramp position, 
bypass door position, angle of attack, angle of sideslip while maintaining 
constant Mach number, airflow level and altitude are required to properly 
assess the data comparisons between flight and wind tunnel data. 

The analysis times for the flight data ranged from .600 to 2.8 seconds 
depending on the desired test condition. Using the acoustic scaling criteria 
of Reference (2), the respective analysis times for the l/6th scale inlet 
wind tunnel model are .100 to .467 seconds. 

Engine Airflow Level 

An engine power setting of intermediate power or afterburning power was 
required for most of the flight maneuvers to achieve or sustain the flight 
maneuver. These two power settings both have maximum corrected engine airflow. 
Therefore, most of the data comparisons for flight test conditions are limited 
to one corrected engine airflow. However, for the wind tunnel test articles, 
the engine power setting was varied or the mass flow plug was translated to 
provide a range of corrected engine airflow levels for both supercritical 
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(without engine only) and subcritical inlet operation for each test condition. 
Therefore, data points for two engine corrected airflow lev 'Is were selected 
from the wind tunnel data to bracket the airflow from the flight test data 
point. This allowed interpolation to match the flight test airflow levels. 
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DATA ACQUISITION AND REDUCTION 

Two similar procedures were used Co acquire and reduce Che time variant 
pressure data, one for the wind tunnel tests and one for the flight tests. 

The primary difference between the two procedures is that the Analog Distor- 
tion Calculator was used in reducing the flight test data but not the wind 
tunnel test data. 

The Analog Distortion Calculator (ADC) was used during the NASA F-15 
inlet/engine compatibility flight test program. Analog values of distortion 
data were computed as a function of time, and electronic flags were placed on 
the wide band FM pressure data tapes corresponding to peak distortion values. 
The analog pressure data were then digitized +300 milliseconds about this peak 
flag., and the digital time variant distortion data were calculated. This 
procedure which was also used during the F-15 inlet development wind tunnel 
program represented a substantial cost savings since it eliminated digitizing 
lengthy data segments (thirty seconds for the F-15 inlet development program). 
However, for this study, the ADC was not used for the wind tunnel data reduc- 
tion procedure since the maximum analysis time was only 2.8 seconds. The 
analog pressure data were digitized, and peak distortion values were calculated 
for the entire analysis time. The data segment to be digitized was randomly 
selected from within the thirty second data record. Approximately 35 million 
pressure values were digitized and used to compute the time variant distortion 
data and the required statistical functions for this study. The data reduc- 
tion paths for the flight test data and wind tunnel test data, the filter 
characteristics, the pressure probe substitution technique, the distortion 
data accuracy and the time history comparisons of analog and digital dis- 
tortion data are described below. 

MCAIR Flight Test Data Processing System 

The system consists of a network of equipment located at MCAIR/Edvards 
and at MCAIR/St. Louis linked together by a high-speed digital transmission 
system utilizing a leased telephone line. The heart of the system is a high 
speed digital computer located at the MCAIR/St. Louis facility. The total 
system is composed of the preprocessor subsystem, digital transmission sub- 
system, high speed digital computer subsystem and the formal data display 
subsystem. A flow diagram depicting the data processing system is shown in 
Figure 9. 

Pr eprocessor Subsystem - The purpose of the preprocessor subsystem is 
to accept the onboard tape recorded signals, perform the necessary processing 
controls, convert the analog signals into a digital form, merge the analog 
and pulse coded modulation (PCM) data, provide time synchronization of the 
data and produce a computer compatible tape. Header information on this 
preprocessed tape provides the necessary instructions to the high speed digi- 
tal computer to allow it to perform all the requested computations and plots 
on the test runs preprocessed. The preprocessor has a mini-computer central 
processing unit (CPU) and a special purpose digital computer. Set-up of the 
preprocessor is accomplished through the use of process control cards that 
include the job definition cards and measurand definition cards. These cards 
are denoted as the Instrumentation Ground Equipment (IGE) cards and are also 
used for pre-flight instrumentation checkout and set-up of the programmable 
airborne data acquisition system. In addition, the preprocessor subsystem 
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uses a data compression algorithm that reduces data transmission time by 
removing redundant data points during transmission. 

Digital Transmission Subsystem - Low- and high-speed digital transmis- 
sion subsystems utilizing voice-grane telephone lines interconnect the MCAIR/ 
St. Louis and the MCAIR/Edwards facilities. The low speed system is used for 
computer file updates and the high-speed system is used as a tape-to-tape 
system managed on both ends by identical controllers. The high speed trans- 
mission system transmits synchronously (constant bit rate) and operates in 
a full duplex mode which allows simultaneous transmission of data and error 
signals. Effective transmission rates of approximately 67C0 bits per second 
are being obtained. The lew-speed transmission systems at the MCAIR/Edwards 
facility consists of a mini-computer interface subsystem configured with a 
keyboard and printer, card reader, and card punch. It is connected to the 
high speed digital computer in St. Louis. Transmission is asynchronous at 
rates up to 14.8 characters per second in a half-duplex mode. 

High Speed Computer Subsystem - This high speed digital computer subsystem 
provides storage space for the information files and the computational capa- 
bility to rapidly process and display the required data in support of the 
MCAIR flight test activities. To handle this task, the system is configured 
with a 512 thousand byte memory, a 400 million byte disk storage, four nine- 
track dual-density (800/1600 BPI) tape units, a line printer (1100 LPM) and 
a combination card reader/punch (1100 C?M read and 300 CPM punch). The sub- 
system is operated in a Virtual Storage (VS-1) environment which is currently 
configured to allow four multiprogrammed partitions of 512 thousand bytes. 

If necessary, any partition size can be increased to 4 million bytes of vir- 
tual memory. This arrangement provides many features of a very large comput- 
ing system without the large scale costs. 

Formal Data Display Subsystems - The Formal Data Display subsystem con- 
sists of line printers and electrostatic plotter/printers. The MCAIR/St. 

Louis facility has an electrostatic plotter/printer interfaced (on-line) 
to the high speed digital computer. This plotter/printer is capable of pro- 
ducing a completely annotated hard copy plot in approximately 4 seconds and 
a tabular printout at approximately 2500 lines per minute. A printer (1100 
lines per minute) is also interfaced to the high speed digital computer. The 
MCAIR/Edwards facility has an electrostatic plotter/printer which operates 
as a "stand alone" system controlled by a mini-computer. In this configura- 
tion, a magnetic tape is the input to the subsystem. 

MCAIR Wind Tunnel Test Data Processing System 

The data processing system consists of a network of equipment located 
at MCAIR/St. Louis. The total system is composed of a data recording sub- 
system, data digitizing subsystem, a high speed digital computer subsystem 
and a formal data display system. 

Data Recording and Playback Subsystem - The high response wind tunnel 
pressure data is recorded on one inch wide band magnetic tape using constant 
bandwidth (CBW) techniques. The subscale data were recorded with six data 
channels per track and the full scale data were recorded with seven data 
channels per track. In addition, two tracks were allocated to record the 
IRIG time code and a reference oscillator signal, and two tracks were left 
blank for use during the digitizing process. 
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Data Digitizing Subsystem - To identify the data segment for digitiza- 
tion, the magnetic tapes are replayed on the tape recorder as shown in Figure 
10a. During playback, the IRIG time code is electronically processed by the 
time code reader and associated electronics and a D.C. level digitizing pulse 
is recorded on one of the blank tracks to identify the data segment selected 
for digitization. 

The tapes were then replayed again as shown in Figure 10b, and the data 
were discriminated (two cracks for each pass), filtered with the desired four 
pole Bessel filters and digitized using sample and hold modules and an 
analog-to-digital converter. Since the data were digitized during multiple 
passes of the tape, it was necessary to time correlate these data. This 
was accomplished by gating the D.C. level digitizing pulse with the IRIG 
time code and then gating the reference oscillator signal to the analog-to- 
digital converter. The digitized data were then buffered into a core memory 
and recorded on a nine track, eight hundred bits per inch digital tape. 

High Speed Digital Computer Subsystem - Using the digital data tape 
generated during the digitizing process as the primary input, the multiple 
pass digitized data segments were merged into a single pass data segment 
which was output to magnetic tape (merged high response data) having a blocked 
format with one data file per data segment, as shown in Figure 10b. The merg- 
ing process was controlled by a high speed digital computer which used a 
80 megabyte disc for working files. The output data format was a 24 bit 
binary integer which corresponded numerically to the 16 bic input data. 

The merged high response pressure data were then read from the magnetic 
tape by a high speed digital computer and calibration scaling was applied to 
obtain high response pressure data in engineering units. These data were 
used to compute the time variant fan distortion data. The time variant pres- 
sure and distortion data were then output to magnetic tape. 

Formal Data Display Subsystem - Using the time variant pressure and 
distortion data magnetic output tape and the appropriate software packages, 
high quality time history traces of distortion data, total pressure contour, 
power spectral density plots, cross spectral density plots, autocorrelation 
plots and cross correlation plots were produced using a high speed 
electrostatic printer/plotter. 

Filter Characteristics 

Low band pass, 4 pole adjustable Bessel filters were used for this study, 
the frequency response characteristics of which are shown in Figure 11. The 
filter cutoff frequency could be adjusted, thus allowing the entire set of 
data used for the study (wind tunnel and flight test) to be reduced using 
this one set (48) of filters, thereby, eliminating inconsistencies in fre- 
quency response characteristics for the model and flight data. High response 
pressure data used to compute time variant distortion data for the F100 engine 
stability audits are filtered at 170 Hertz cutoff frequency as specified 
by P&WA which corresponds to design rotor speed. 
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Pressure Probe Substitution Technique 

During the computation of digital distortion and generation of engine 
face pressure profile plots, it was sometimes necessary to substitute for 
inoperative probes on the rake to otain values at all 48 probe locations. 

A probe was considered to be out of tolerance when the total pressure 
recovery for that probe location was greater than 1.1 and less than 0.2. The 
output from each probe of the rake was examined prior to data reduction, and 
probes were deleted if the probe output was suspect. In addition, data 
reduction was halted if more than six probes were inoperative on two adjacent 
legs or more than twelve probes were inoperative on the forty eight probe 
rake. 


A set of probe substitution criteria which were developed during the 
F-15 inlet envelopment program were used for the study are are presented in 
Figure 12. If these criteria failed to determine a forty-eight probe set, 
the adjacent most radial and circumferential three probes were averaged and 
used to obtain a probe location valu Probes were obtained by proceeding 
in a radial direction first and in a rcumferential direction second. 

Distortion Data Accuracy 

An error assessment was made in Reference (3) for the F-15 flight test/ 
full scale inlet wind tunnel model engine face rake configuration. Typical 
error in the fan distortion factor as a function of pressure error is shown 
in Figure 13. 


Analog and Digital Distortion Data Comparisons 

The Analog Distortion Calculator was used during the NASA F-15 inlet/ 
engine flight test program to flag peak values, therefore identifying the 
data time segment to be digitized. It is important that the ADC calculate 
peak time variant distortion accurately so that the correct time segment is 
selected for digitization. To illustrate the accuracy of the ADC, time his- 
tory traces of analog distortion data and of digital distortion data for the 
flight test program are compared in Figure 14 and additional comparisons 
are made in Appendix A. As can be seen from the figures, the ADC computes 
peak distortion values to within +5 percent of the digitally computed peak 
distortion values. 
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DATA BASE SELECTION 


The data base for this stu^y was established by assembling data from 
four distinct inlet test series: (1) the NASA inlet/engine compatibility 

flight test program, (2) the MCAIR full scale inlet development test without 
engine (FSCP) , (3) the MCAIR full scale irlet development test with engine 
(FSE), and (4) the MCAIR l/6th scale inlet development test. To meet the 
objectives of the study, test conditions were selected which encompassed 
wide ranges of Mach number, angle of attack, angle of sideslip and corrected 
engine airflow and had high levels of peak distortion. 

The inlet configuration, as previously described, was identical for all 
the wind tunnel models and the flight test vehicle. The programmed inlet 
ramp and bypass door scheuules were also identical. The inlet ramp angles, 
as defined in Figure 2, were identical to within +0.5 degrees for each 
selected data point and test article for a given test condition. The Dypass 
door areas were either fully open or fully closed for each test condition 
with the exception of three data points which had partially opened bypass 
positions. The test conditions for each data point, model scale, Mach 
number, angle of attack, angle of sideslip, inlet cowl rotation angle, 
inlet third ramp angle, bypass door position, percent corrected engine air- 
flow, Reynolds number, analysis time and part-point number for this study 
arc summarized in Figure 15. 
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DISTORTION DESCRIPTOR 


One of the objectives of this study was to evaluate the capability of 
the Pratt & Whitney stability audit procedure to predict engine stalls. 

Since current state of the art distortion descriptors are dependent upon the 
engine configuration, evaluation of the Pratt & Whitney stability audit 
procedure necessitates using the Pratt & Whitney distortion descriptor, Ka 2 . 
since it is the distortion descriptor for the F-15 poverplant, the F100 
engine. 

During the development phase of the F-15 inlet/engine propulsion system, 
the distortion descriptor was repeatedly improved as the experimental data 
base increased to better define the relationship between the distortion 
descriptor and the engine stall margin. As a result of the repeated improve- 
ment of tne distortion descriptor, the wind tunnel distortion data and flight 
test distortion data were not comparable on a direct basis. To make valid 
comparisons of distortion levels, the distortion data for all the test 
conditions of Figure 15 were computed using the current distortion descriptor 
per Reference (4). The oat nratical formulae of the distortion descriptor 
used to reduce the data for this study and previous F-15 development data are 
presented in Figure 16. The parameters that have been changed throughout the 
development phase of the F100 engine are the reference radial profile and 
the radial distortion weighting factor. A sample calculation of the distor- 
tion descriptor is shown in Appendix 3, and the reference radial profiles 
and radial distortion weighting factors for the F100 (2 7/8) engine are also 
presented. 

Distortion descriptor values are derived to have a linear relationship 
with remaining engine stall margin, and the distortion limit is defined as 
the distortion level for which no remaining stall margin exists (engine stall). 
Therefore, the significance of a distortion level is its relative magnitude 
to the distortion limit ievel for which an engine stal^ may occur. The 
levels of distortion values for which a stall may occur for the F100 (2 7/8) 
engine are shown in Figure 17. This curve is a screening curve which is con- 
servative relative to predicting engine stalls since it uses the lowest 
level of distortion for a given airflow that could possibly stall an engine 
for any given Mach number and altitude. It is termed a screenirg curve since 
values from this curve are not to be used as distortion limits for engine 
stall; instead when an inlet distortion level for a particular airflow ex- 
ceeds the value from this curve, additional analysis should be conducted to 
determine if there is any remaining engine stall margin for the particular 
test condition. It should be noted from Figure 17 that the variation in 
engine tolerance to distortion level from high to low airflow is significant 
since the respective ranges of distortion level from the screening curve 
range from 0.75 (high airflow) to 16.75 (low airflow). 

The range in distortion level from high to low airflow is primarily a 
result of significant changes in the radial distortion weighting factor, b. 

This factor is shown in Figure 18 to range from 0.93 (high airflow to 29.29 
(low airflow) as a function of airflow. Since fan distortion is the sum of 
circumferential distortion and radial distortion multiplied by the radial 
distortion weighting factor (Ra^ ■ K + b • K r ) , and since the levels of 
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circumferential and radia distortion generally have values between 0.25 to 
1.0 for the engine airflows of interest in this study, radial weighting fac- 
tors corresponding to low engine airflows (b * 25) can easily result in fan 
distortion values of 16 or greater. The radial distortion weighting factor 
implies that at reduced engine airflows, the fan is more sensitive to radial 
distortion than to c ire inafe rent ial distortion. As a result of discontinui- 
ties in the radial distortion weighting factor. Figure 18, at approximately 82 
percent corrected engine airflow, an interpolation procedure was required 
for distortion data comparisons at constant corrected engine airflows. This 
is discussed in the Interpolation Procedure section. 

The irregularity in the radial distortion weighting factor is directly 
related to the subtle aspects of the P&WA F100 stability audit system. The 
decreasing sensitivity of the F100 fan to circumferential distortion and the 
increasing sensitivity to radial distortion, as corrected engine airflow is 
reduced, results in the irregularity at approximately 82 percent corrected 
engine airflow as shown in Figure 18. 


17 



TIJIE VARIANT DATA QUALITY 


The aata comparisons necessary to satisfy the primary objectives of 
this study require that the samples of time variant distortion data be 
stationary random data. This implies that the statistical properties of 
each data sample are constant with time and are repeatable. With constant 
statistical properties and identical test conditions, reduced data from the 
various inlet test programs can be meaningfully investigated for the effects 
of Reynolds /scale effects, frequency content and engine presence. Detailed 
investigation into the data classification, description and characterization 
were made and samples of wind tunnel and flight data were evaluated in detail 
for stationary random data verification. In additior, comparisons of sta- 
tionary random data and periodic data were made for two samples of flight data. 

The test conditions selected from the complete data matrix. Figure 15, for 

evaluation of stationary random characteristics are shewn in Figure 19. 

Data Classification 

Data from observed phenomena are generally categorized as either deter- 
ministic or random, and the various classifications of data for both deter- 
ministic and random data are shown in the schematic of Figure 20. Determin- 

istic data can be categorized as being either periodic or non-periodic. 

Periodic data can be further categorized as being either sinusoidal or non- 
sinusoidal (complex periodic). Non-periodic data can be further categorized 
as being either "almost periodic" or transient. In addition, combinations 
of these deterministic forms may also occur. 

Data representing random physical phenomenon cannot be described by an 
explicit mathematical relationship siac^ each observation of the phenomenon 
is unique. Any given observation will represent only one of many possible 
results that could have occurred. A single time history representing a ran- 
dom phenomenon is called a sample record, and the collection of all possible 
sample records which the random phenomenon could have produced is called a 
random process. Therefore, a sample record can be thought of as one physical 
occurrence of the random process. Random processes can be categorized as 
either non-s tat ionary or stationary. Non-s tat ionary random processes can 
be further categorized into specific types of non— scat ionary phenomenon. 

More detailed discussions of non-s tat ionary data can be found in Reference 
(6). Stationary random processes can be further categorized as either 
ergodic or nonergcdic. 


Data Description 

Examples of each data classification are as follows: 
o Deterministic Data 

- Sinusoidal Periodic - Described mathematically by the periodic function: 
x(t) * X sin (Z^fo + 0) 

- Complex Periodic - With few exceptions can be described mathematically 
by the Fourier series: 
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x(t) ■ ■— + z [a Q cos (2rnf^t) + sin (2xnf^t)j 
n*l 

- Almost Periodic - Has an infinitely long fundamental period and does 
not satisfy: x(t) - [x(t) + NT p ] relationships. 

Example: x(t) * sin (2t + 6^) + X, sin (3t + 6^) 

+ X, sin (^50t +3,) 


o Random Data 

- Mon-Stationary - includes all random processes which do not meet the 

requirements for stationarity. Example: x(t) * X(t) • 

A(t) where X(t) is random and A(t) is deterministic. 

- Weakly Stationary - Mean Value (first moment) invariant with time 

- Autocorrelation (joint moment) invariant with time 

- Strongly Stationary - All possible moments and joint moments invar- 

iant with time. 

- Ergodic Stationary - Mean value (first moment) invariant with time 

for all sample records 

- Autocorrelation (joint moment) invariant with 
time for all sample records 

For deterministic data in general, explicit mathematical functions will 
describe the dependent variable as a function of time for a particular cate- 
gory of data. In contrast, there are no explicit mathematical functions 
that will describe the random data observation since it is unique. However, 
time averages of stationary random data will be constant over the sample 
record of interest, and if the data is ergcdic stationary random, the time 
averages will be constant for all sample records for a given ensemble of 
data. In actual practice, if random data exhibit stationary physical phenom- 
ena, the data are generally said to be ergodic and random stationarity can 
be measured properly from a single sample record rather than an ensemble of 
simple records. 


Data Characterization 

Since no explicit mathematical functions exist for stationary random 
data, four main types of statistical functions are used to characterize the 
basic properties of stationary random data. These statistical functions are 
(a) mean square values, (b) probability density functions, (c) autocorrelation 
functions, and (d) power spectral density functions. 

The mean square value furnishes a description of the intensity of the 
data. The probability density function (PDF) furnishes information about the 
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data in the amplitude domain. The autocorrelation function and the power 
spectral density function furnish information about the data in the time 
domain and the frequency domain respectively. For stationary random data, 
the power spectral density function furnishes no new information relative to 
the autocorrelation function since the autocorrelation function and the power 
spectral density function are transform pairs of one another. 

For most stationary random data, it is unnecessary to derive all of the 
statistical functions discussed above. Csing four basic assumptions, which 
are listed below, stationary random data can be verified using the time 
varying mean and time varying mean square values. These statistical param- 
eters must have a constant value for the sample record. It should be noted 
that the mean square of the fluctuating pressure component is equal to the 
variance of the instantaneous pressure (steady state pressure component plus 
the fluctuating pressure component) and is referred to as the instantaneous 
pressure variance throughout the remainder of the report. 

T T 

Mean Square of . . . ? ~ 

Fluctuating Pres- * ^ - f [ P^U)]" dt = y / [P Ci (t) - p t ss l dt 
sure Component 

o o 

= Variance of Instantaneous 
Pressure 

Theoretical proof of stationarity requires that all statistical properties 
of the random process be invariant with time. Verifications of this type 
are clearly not feasible since there are an infinite number of statistics that 
would have to be computed for an ensemble rather than a single sample 
record. However, four basic assumptions can be made, thereby allowing prac- 
tical tests for stationarity. 

The first assumption is chat the statistical properties of a short time 
interval of data will not vary significantly from a short time interval of 
the next sample. This simply means that self stationarity of individual 
records can be accepted as stationarity of the random process from which the 
sample records were taken. 

The second assumption is that verification of weak stationarity is suffi- 
cient. If this assumption is accepted, stationarity can be verified using 
only the mean value and the autocorrelation function of the sample record. 

Two reasons make this assumption acceptable. Power spectra and autocorre- 
lation analysis do not require more than weak stationarity to be valid and 
random data will generally be strongly stationary if weakly stationary. For 
the case of data with a normal probability density function, weak stationary 
data are strong stationary data since all higher order statistical properties 
of normal data are determined by the mean value and autocorrelation function. 

The third assumption is that the sample record to be investigated is long 
compared to Che random fluctuations of the data time history. In effect, the 
sample must be lengthy enough to permit long term trends to be differentiated 
from the time history fluctuations. 
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A fourth important assumption is that if the variance of the sample data 
record is stationary then che autocorrelation function is also stationary. 

This assunption which is somewhat less dependable than the first three 
assumptions simplifies the testing procedures and is usually valid since 
it is highly unlikely that nonstationary data will have a time varying auto- 
correlation function for any time displacement T without the value at t * 0 
varying. Since the variance equals the autocorrelation function at t * 0, 
the variance racher than the entire autocorrelation function is used. 

With these assumptions, the stationarity of random data can be verified 
using only the mean and the variance of the sample data record. These assump- 
tions are verified in the following section by examining power spectral 
density functions, autocorrelation functions, probability density functions, 
mean values and variance values for selected flight test and wind tunnel 
data points. 


Verification of Random Stationary Data 

Statistical data analyses were made to: (1) verify the use of the pres- 

sure variance and the mean pressure values to determine if data are random 
stationary, (2) verify the flight test data to be stationary random for 
relatively short analysis times and (3) illustrate the differences between 
stationary random data and periodic (deterministic) data. 

Oue of the most extreme flight test conditions, Mach number 1.6 at angle 
of attack of -4 degrees and sideslip angle of 0 degrees, was selected to 
verify the assumptions for stationarity of random data. This test condition 
required a transient flight maneuver, shown in Figure 8, to obtain the desired 
flight attitude. The test condition was obtained by a pull-push maneuver and 
the aircraft was on condition for approximately 0.6 seconds as noted in Figure 
15. If this data point, with the underlying assumptions for stationarity 
exhibits the properties of stationary data, che assumptions for stationarity 
should also be valid for the wind tunnel data and one "g" flight conditions. 

The power spectral density function (PSD) of Figure 21(a), generated 
using data for the flight test condition from Figure 8, exhibits the same 
frequency characteristics of wide band random noise typical of random sta- 
tionary data, per Reference (6). In contrast, the PSD function of Figure 
21(b) shews periodic data with discrete frequency spectra of approximately 
9, 18, 27, 39, 45 and 54 Hertz. This contrast of the PSD functions empha- 
sizes the pressure data randomness (no discrete frequency spectra) of the 
selected flight test condition. 

The autocorrelation functions for the two flight test conditions are 
presented in Figure 22. The autocorrelation function of Figure 22(a) shows 
no correlation for the inlet rake pressure probe L6R3 at the engine face. 

This is typical of random data, and very rapidly approaches a constant value 
as a function of lag time which also is typical of stationary data, there- 
fore the data again exhibits the properties of stationary random data. For 
the periodic data of Figure 22(b), almost perfect correlation is obtained 
for the same pressure probe location as used for the random data analysis. 

The elapsed times between the peak correlations for the pressure probe of 
Figure 22(b) are 0.111 seconds which correspond to the fundamental frequency 
of 9 Hertz obtained from the PSD function. As expected, the autocorrelation 
function does not approach a constant value. 


21 



The probability density functions for the random and periodic data sets 
for the same pressure probe used in the comparisons of the PSD and autocorre- 
lation functions are plotted in Figure 23. Even though only a limited number 
of pressure data are available, the random data set are seen to approach a 
normal distribution and the periodic data set are seen to approach a sine 
vave distribution as expected. Since the random data is approximately nor- 
mally distributed, the higher order statistical properties will be stationary 
since the autocorrelation function (variance when using the fourth assumption 
as described above) and the mean are stationary. 

The time history values of freestream total pressure, instantaneous probe 
pressure L6R3 and the associated mean and variance values are presented in 
Figure 24 for both the random and periodic flight test data sets. The 
time history mean and variance values can be seea to be stationary for the 
random data of Figure 24(a), whereas the time history mean and variance 
values of Figure 24(b) are periodic as expected for the periodic data. 

The statistical functions used to prove random stationarity for the flight 
data have also been used to verify random stationarity for the wind tunnel 
data. These statistical functions are presented in Figure 25 through 
Figure 30 for data collected using the full scale model with engine, the 
full scale coldpipe model (without engine) and the l/6th scale inlet model. 
High frequency spikes occurred in some of the wind tunnel PSD plots as shown 
in Figures 25a, 27a and 29a. The cause of these spikes is unknown, however, 
the frequency at which they occurred is significantly above the frequency of 
interest for this study. The normalized autocorrelation function of Figure 
27 shows a bias or offset of approximately -0.2. The cause of the offset is 
unknown, but the level of correlation is insignificant. These two discre- 
pancies have no impact on this study. The time history plots of instantaneous 
pressure, of mean pressure and of the pressure variance for selected probes 
from each of the test conditions listed in Figure 19 are presented in Appen- 
dix C. 

In summary, representative data from the complete data matrix have been 
shown to be random stationary. In addition, assumptions have been verified 
that allow random stationary data to be demonstrated using only the time 
averaged mean and variance of the pressure values. 
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INTERPOLATION PROCEDURE 


One of the primary considerations in establishing the study base was 
identifying test conditions with comparable engine airflows for each test 
article. Since actual inflight engine airflows could not be predicted for 
every combination of Mach number and altitude combination during the develop- 
ment phase of the inlet, the data collected during the wind tunnel testing 
and the data collected during flight testing generally have different engine 
airflow values even though the geometrical configurations are identical. 
Significant differences in distortion level can occur not only because of the 
differences in the total pressure contours at the engine face, but also as a 
result of significant changes in the engine related weighting factors which 
are a function of engine airflow. 

As an example, the effect of filter cutoff frequency on fan distortion 
and turbulence is illustrated in Figure 31 where the airflow for each data 
point of the test condition is different. These data are not directly com- 
parable since each curve is for a different airflow. Therefore, to illustrate 
the effect of engine airflow, the data with filter cutoff frequency of 170 
Hertz are plotted in Figure 32 as a function of corrected engine airflow for 
each of the three test articles. Since only one airflow level (78.9 percent) 
is available for the flight data, the full scale cold pipe data and the full 
scale with engine data have to be interpolated and extrapolated respectively 
to obtain comparable distortion and turbulence values. Note that the full 
scale cold pipe fan distortion levels at 78.9 percent engine airflow are 
significantly higher than those of the other two test articles. 

In general, the circumferential and radial components of fan distortion 
are well behaved functions of airflow. Since fan distortion is not only 
dependent on these descriptors but also on the radial distortion weighting 
factor which is not a well behaved function of airflow (as discussed above), 
fan distortion is not expected tc be a well behaved function of airflow. 
Therefore, interpolation of circumferential and radial distortion data is 
technically feasible whereas interpolation of fan distortion is questionable. 

A technique involving linear interpolation/extrapolation of circumfer- 
ential and radial distortion data was developed for this study. The inter- 
polated values are used together with the appropriate radial distortion 
weighting factor as a function of engine airflow to define fan distortion 
values at constant engine airflow. The solid symbols of Figure 33 were 
computed using this technique. In the engine airflow range of 82 percent, the 
interpolated values of fan distortion strongly reflect the characteristic 
variation of the radial distortion weighting factor with airflow, as might be 
expected. However, comparison of the differences between the distortion 
levels (j4).4) at 78.9 percent for the three test articles and of the screening 
curve. Figure 17, distortion level (8.4) illustrates that the differences 
are relatively insignificant. 

The results of interpolating the fan distortion and turbulence at ail 
frequencies for 78.9 percent engine airflow is shown in Figure 34. For a 
constant filter cutoff frequency, increasing turbulence correlates with in- 
creasing fan distortion, and for a particular test article increasing the 
filter cutoff frequency increases both fan distortion and turbulence. 
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Interpolated/extrapolated fan distortion values with respect to airflow 
have been computed using this technique throughout the study to make direct 
data comparisons when comparisons at identical airflows are not available. 
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EXCISE PRESENCE EFFECT ON DISTORTION LEVELS 


The data for this study are unique since there are limited sets of 
data where the same inlet has been tested both with and without an engine. 

The effect of the engine installation was evaluated for the test conditions 
shown in Figure 35. Sixteen data points for supersonic flight conditions 
were used to establish the five test conditions for evaluation. 

The inlet models used to generate the data are the full scale with 
engine (FSE) and the full scale cold pipe (FSCP) which in effect is the 
same inlet model tested with and without an engine. Therefore, the inlet 
geometrical configurations are identical for these test conditions with 
the exception of bypass door areas which differ for two test conditions as 
noted in Figure 35. 

However, the bypass area differences were only significant for the 
Mach 2.2, a = -2° test condition, therefore only this condition was elimin- 
ated from the comparisons to determine engine presence effect. In deter- 
mining the effect of engine presence on distortion level, the following six 
items were investigated: (1) steady state and instantaneous pressure con- 

tours, (2) power spectral density functions, (3) derived turbulence level, 

(4) total pressure recovery, (5) fan distortion descriptor, and (6) mea- 
sured turbulence level. 

Contour Plots and Power Spectral Density Plots 

The comparisons of steady state total pressure contours, instantan- 
eous total pressure contours and power spectral density (PSD) plots with 
and without engine are shown in Figures 36 through 45. In general, the 
steady state total and instantaneous pressure contours compare favorably. 

In particular, the steady state contours of Figure 40 are almost identical. 
The engine airflow differences between these contours is 1.0 percent. The 
small differences in the steady state contours of the other comparisons may 
be due to small differences in the airflow levels. 

Although the patterns for with and without engine are similar for both 
steady state and instantaneous pressure contours, one qualitative difference 
was observed. In Figure 40, the "islands" in the upper right quadrant for 
both the steady state and the instantaneous contours tend to be slightly 
rotated counter clockwise for the with engine contours, the same rotational 
direction as the engine. 

The steady state pressure contours are labeled with the static pressure 
ratios and their locations in Figures 36 through 40. A compilation of these 
static pressure ratios is presented in Figure 36 for increasing engine air- 
flow for each test condition. In addition, the average hub and average tip 
static pressure ratios have been computed for each data point and are plotted 
in Figure 47. For these five test conditions, the data indicates that the 
engine presence decreases the average tip static pressure ratio relative 
to the average hub static pressure ratio. These relative changes in hub and 
tip static pressure ratios reflect relative changes in hub and tip duct velo- 
cities that are conducive to improved inlec engine compatibility for these 
test conditions. 
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The power spectral density (PSD) plots, with and without engine are 
presented in Figures 41 through 45. Additional PSD plots for each of these 
data points are presented in Appendix D, and the corresponding auto-correlation 
functions are presented in Appendix E. The PSD functions are well behaved 
for all the data points and test conditions evaluated. 

Turbulence, derived from these PSD plots by integration, is presented 
as a function of filter cutoff frequency in Figures 48 through 52 for the 
respective PSD plots of this section. The characteristic function of tur- 
bulen. e as a function of frequency is consistent and well behaved for these 
pressure probes. Turbulence values for every PSD plot of this study were 
derived and presented as a function of frequency in Appendix F. No apparent 
engine effect was observed for the PSD plots or for their respective integrat- 
ed functions. 


Pressure Recovery, Distortion and Turbulence 

Total pressure recovery, peak time variant distortion and turbulence 
values with and without the engine present are compared in Figures 53 and 
54. In these figures, favorable engine effects are indicated for total 
pressure recovery values above and time variant distortion values below 
the line of perfect agreement. For each of these three inlet parameters, 
the data indicated that the engine has a small but favorable effect. 

In summary, a qualitative assessment of the total pressure contours 
and power spectral deusiLy plot comparisons indicates that the "with" and 
"without engine" inlet configuration result in no significant differences. 
However, a quantitative assessment of total pressure recovery, peak time 
variant distortion and turbulence level comparisons indicate a small but 
favorable effect of the engine installation. In addition, the data indicate 
that the engine presence decreases the average tip static pressure ratios 
relative to the average hub static pressure ratios. 
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REYNOLDS NUMBER/SCALE EFFECTS ON DISTORTION 
PRESSURE RECOVERY AND TURBULENCE 


To evaluate the effect of Reynolds number/ scale on peak time variant 
distortion, steady state distortion, total pressure recovery, spatial dis- 
tortion and turbulence, 62 data points were utilized from a total of 82 for 
the study. The respective test conditions for the 62 data points are shown 
in Figure 55. The Reynolds number/scale evaluation was conducted using 
time variant data filtered at 170 Hertz for full scale data and at 1040 Hertz 
for l/6th scale data. The Reynolds number/scale values are computed by 
multiplying the Reynolds number per foot values obtained from the test 
conditions by one foot (unit characteristic length) for full scale data and 
by l/6th foot for the l/6th scale data. 

In general, the range of Reynolds number/scale values for each of the 
test articles are shown in the following table. An F-15 flight envelope 
as a function of Reynolds number/scale and Mach number is shown in Figure 56. 


Inlet Model 

Reynolds Number/Scale 

l/6th scale 

0.2 x 10 6 - 0.4 x 10° 

Full Scale Cold Pipe (FSCP) 

1.2 x 10 6 - 1.5 x 10 6 

Full Scale With Engine (FSE) 

1.2 x 10 6 - 3.6 x 10 6 

Flight 

0.8 x 10 6 - 3.6 x 10 6 


Plots cf total pressure recovery, turbulence levels, peak time variant 
distortion values and time variant spatial distortion (at peak time variant 
distortion) as a function of Reynolds number/scale are presented in Figures 
57 through 82. In addition, steady state fan distortion data and steady 
state spatial distortion data are shown as solid symbols and exhibit the 
same trends as uhe corresponding time variant data. The majority of the data 
show the expected trends of increasing total pressure recovery and decreas- 
ing peak time variant distortion and turbulence with increasing Reynolds 
number- scale. Time variant spatial distortion was evaluated at the setae time 
frame at which the peak time variant distortion level occurred and shows no 
consistent trend with increasing Reynolds number/scale. 

The data trend of increasing pressure recovery and decreasing fan dis- 
tortion and turbulence as a function of increasing Reynolds number/scale does 
have some inconsistencies, examples of which are shown in Figures 58 and 62. 
The discrepancies occur in comparing l/6th scale inlet data at low Reynolds 
number/scale values (.2 - .45 x 106) to full scale inlet data (either flight, 
FSE or FSCP) at low Reyolds number/scale values (approximately 1.5 x 1C)6) . 
These discrepancies do not occur when comparing l/6th scale inlet data at 
low Reynolds number/scale to full scale data at large Reynolds number/scale 
values (3.5 x 10^). The discrepancies noted represent only a small percent- 
age of the data evaluated for the effect of Reynolds number/scale and may 
be due to several possible factors. 
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Peak time variant distortion is a statistical parameter since it is 
computed from random stationary time variant pressure data. The peak value 
for a given time segment may not repeat itself exactly, but may fall into a 
range of values for the specified time segment. Therefore as the range of 
Reynolds number/scale values is reuuced, the variation in the peak distortion 
values for a given time segment may overshadow the change in peak time 
variant distortion for small changes in Reynolds number/scale values. Note 
that the discrepancies did not occur for large differences in Reynolds number/ 
scale values. In addition, the discrepancies in the Reynolds number/scale 
data trend occurred for comparisons of l/6th scale inlet data and full scale 
inlet data. While these two sets of data were acquired with identical rake 
configurations, different dynamic probes and different data acquisition 
systems were used. 

In summary, the data exhibits a definite trend as a function of Reynolds 
number/scale, i.e., increasing total pressure recovery, decreasing peak time 
variant distortion values and decreasing turbulence levels as a function of 
increasing Reynolds number/scale. The time variant spatial distortion value 
take., at peak time variant fan distortion did not correlate with the peak 
time variant fan distortion value. 



FREQUENCY CONTENT EFFECT ON DISTORTION AND TURBULENCE 

The test conditions and the thirty two associated data points used to 
determine the effect of frequency content on peak time variant fan distor- 
tion, turbulence level and time variant spatial distortion values (at peak 
time variant fan distortion) fot various values of Reynolds number/scale are 
shown in Figure 83. For each test condition, the data were evaluated at the 
engine corrected airflows listed in Figure 83. 

To relate frequency content to peak time variant distortion, the analog 
pressure data were filtered at selected filter cutoff frequencies, and peak 
time variant fan distortion, turbulence and spatial distortion were calcu- 
lated and plotted as a function of the filter cutoff frequency. The filters 
used for this study were the adjustable four pole Bessel filters described 
in Figure 11. The full scale filter cutoff frequencies were selected to be 
both greater than and less than the P&WA specified engine filtering frequency 
of 170 hertz. The l/6th scale filter cutoff frequencies are scaled to match 
the full scale filter cutoff frequencies using the frequency scaling criteria 
of Reference (2). The filter cutoff frequencies used for full scale pressure 
data were 45, 100, 170 and 500 hertz, and the filter cutoff frequencies used 
for the l/6th scale pressure data were 275, 615, 1040 and 3070 hertz. The 
frequency content effect on fan distortion, turbulence and spatial distortion 
for the l/6th scale test is plotted as a function of the equivalent full scale 
filter frequencies to simplify the comparisons between the l/6th scale and 
'ull scale data. To accurately represent these data digitally, the analog 
pressure signals were sampled at approximately five times (or greater) the 
appropriate filter cutoff frequency. The sampling rates varied from 270 
samples per second to 21,000 samples per second, and these rates are listed 
in Figure 11 for each of the respective filter cutoff frequencies. 

^ Filter Cutoff Frequency Effect 

Increasi^l filter cutoff frequency is seen to -L.icrease peak time vari- 
ant distortion as shown in Figure 34. This effect of filter cutoff frequency 
on peak time variant fan distortion and turbulence is expected if a physical 
understanding of the frequency content is recognized. A power SDectral 
densitjs&plot is shown in Figure 84a, with dependent and independent parameters 
of AP^jhz and frequency respectively. The area beneath the curve from zero 
frequency to a selected "cutoff" frequency is proportional to turbulence. 

Since turbulence is proportional to the area beneath the PSD curve, the large 1 " 
the cutoff frequency for a particular pressure probe tha large the turbulence 
as shown in Figure 84b. This same logic follows for peak time variant fan 
distortion values. Since tine variant distortion is calculated using time 
variant pressure joiner eases in time variant pressure variance (increases in 
turbulence) increase the probability of larger peak time variant distortion 
values with increasing cutoff frequency. 

Peak time variant fan distortion increases as filter cutoff frequercy 
is increased, but the peaks do not necessarily occur at the same segment time. 
This effect can be seen in Figure 85 where time history traces of time variant 
6 *' .ortion for cutoff frequencies of 45, 100, 170 a d 500 Hertz are shown for 
data point 43. Time history traces of time variant distortion are shown in 
Appendix G for all of the data points in this section. The time history data 
traces include 100 milliseconds of l/6th scale (600 milliseconds for full- 
scale data) time variant distortion data and the peak distortion level. 
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Peak time variant distortion, turbulence and spatial distortion data 
are presented as a function of full scale filter cutoff frequency in Figures 
86 through 93. The curves for turbulence pass through zero turbulence at 
zero cutoff frequency. This simply implies that at zero frequency response 
the tine variant pressure (wit*, no variance) is equal to the steady ^tate 
pressure. Similarly, the time variant distortion value at zero frequency 
response is equal to the steady state fan distortion. The steady state fan 
distortion levels are denoted by solid symbols in the plots. 

For all test conditions, increasing filter cutoff frequency increased 
both the levels of peak time variant fan distortion and turbulence. For 
all but one test condition, the test article with the highest turbulence also 
has the highest distortion level for a given filter cutoff frequency. The 
time variant spatial distortion data at peak time variant fan distortion did 
not indicate any correlating properties with peak fan distortion. The plotted 
values of spatial distortion appear to occur randomly about the steady state 
level . 


Combined Effect of Filter Cutoff Frequency and Reynolds Xumber/Scale 

The effect of Reynolds number/scale on peak time variant fan distortion, 
turbulence level and spatial distortion is shown in Figures 94 through 101 
for constant levels of filter cutoff frecuency. The trend of decreasing 
turbulence level and decreasing peak time variant fan distortion for increas- 
ing Reynolds number/ scale is similar to that shown in the section for Reynolds 
number/scale effects. This trend in the data appears to be consistent for 
ell levels of filter cutoff frequency. Since much of the same data is 
utilized here as in the Reynolds number/scale analysis, some of the same dis- 
crepancies occiu. Since they were discussed in the previous section, they 
are not reiterated here. 

In summary, turbulence level and peak time variant fan distortion in- 
crease as filter cutoff frequency is increased. This data trend is main- 
tained when evaluating turbulence level and peak time variant distortion as 
a function of Reynolds number/scale for constant filter cutoff frequency. 
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PEAK TIME VARIANT FAN DISTORTION PREDICTIONS 


The primary objective of this analysis was to compare the predicted most 
probable peak distortion values as computed by Melick's procedure to the mea- 
sured peak time variant distortion values. Melick's procedure was developed 
under NASA Contract NAS2-6901 and is described in detail in Reference (1) . 

The necessary inputs to Melick's procedure, an evaluation of the dis- 
tortion prediction capability and an assessment of the hypothesis used in 
developing the procedure are described below. 

Inputs to Melick's Procedure 

The necessary inputs to Melick's procedure are: 

o Unfiltered turbulence for each of the dynamic probes to be used 
in the prediction (each data point was evaluated using unfiltered 
turbulence levels for all 48 dynamic probes) , 

o Ratio of filtered mean square pressure variation to unfiltered 
mean square pressure variation for each of the dynamic probes 
to be used in the prediction (each data point was evaluated using 
mean square pressure ratios for all 48 dynamic probes) , 

o Filter cutoff frequency for the filtered mean square pressure varia- 
tion, 

o Airflow velocity at the engine face, normalized dynamic pressure, 
(q/P^) at the engine face and engine cutoff frequency, (f c / 0 ) eng, 

o Reference radial profile, radial distortion weighting factor, and 
analysis time, 

o Steady state total pressure array at the engine face and the 
average engine face static pressure ratioed to rreestream total 
pressure. 

Measured and Predicted Peak Instantaneous Distortion Value Comparisons 

The data points used to evaluate Melick's procedure are shown in 
Figure 102. This data set is composed of two sets of data, one with nominal 
average turbulence with no engine stalls and one with a high levei of average 
turbulence such that an engine stall occurred for each data point. The com- 
parison of measured peak time variant distortion values and the predicted 
most probable peak distortion values for the data with nominal turbulence 
levels and no engine stalls is shown in Figure 103. The average error for 
th^ predicted data is 11.3 percent. The comparisons of measured peak time 
variant distortion values and the predicted most probable peak distortion 
values for the data with high turbulence levels and engine stalls is shown 
in Figure 104. The average error for the predicted data is 20 percent with 
the maximum error approaching 40 percent. Note that the predicted distortion 
levels are lower than the measured distortion levels for all test conditions 
for the distortion data with high turbulence levels. 
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Vortex Diameter and Cross Correlation Coefficient Comparisons 

Melick hypothesized in Reference (1) that inlet turbulence results from 
a random distribution of discrete vortices being convected downstream by the 
mean flow, and in Reference (7), that the cross correlation coefficient (at 
zero lag time) between two dynamic probes is a function of the probe spacing 
normalized by the vortex radius for these set of random vortices. To in- 
vestigate these hypotheses for every pressure probe at the engine face would 
incur an exorbitant data reduction cost. However, Melick indicates in Reference 
(1) how the mean vortex position can be found, therefore allowing judicious 
selection of probes to be cross correlated. The procedure for finding the 
position of the vortex, an example of its use, and an investigation of Melick' s 
hypothesis regarding the presence of vortices are discussed below. 

Using the steady state total pressure data at the engine face, the rake 
legs having the minimum and maximum total pressure recovery are determined. 

If these two legs divide the engine face into two unequal circular segments, 
the vortex should be found in the smaller circular segment midway between the 
two legs. If the two legs with minimum and maximum total pressure recovery 
divide the engine face into two semi-circles, che vortex should be found mid- 
way between the two legs in the semi-circle with lowest overall pressure 
recovery . 

An example of how the probes were selected for computing cross correla- 
tion functions uses the steady state total pressure contour shown in Figure 
105. The rake legs with minimum and maximum total pressure recovery bracket 
the circumferential section of interest and are noted. The probes are 
selected on the rake leg approximately midway within the designated circum- 
ferential sector. The circumferential section and rake leg location .is noted 
in Figure 105. Cross correlation functions for the selected probes are pre- 
sented in Figures 106, 107 and 108, and the peak correlation coefficients 
between the probes are noted in each figure. Cross correlation functions 
for each of the data points evaluated (without engine stalls) a^e presented 
in Appendix H and the respective cross spectral density plots are presented 
in Appendix I. 

Using Melick’s hypotheses of References (1) and (7) as previously dis- 
cussed, the values of the computed vortex radius using Melick' s analysis 
and the Cross Correlation Coefficients (at zero lag time) derived fro- 
measured pressure data are presented in Figure 109 for selected probes lo- 
cated in the region of the hypothetical random vortices. Using Melick' s 
analysis. Reference (7), the cross correlation coefficients were derived that 
correspond to the computed vortex radius and probe spacing for each of the test 
conditions of Figure 109. The derived cross correlation coefficients ana che 
measured cross correlation coefficients for each set of probes have been 
tested for similarity using the Chi-Square Test. Acceptance of the null 
hypothesis (H Q : derived cross correlation coefficient = measured cross cor- 

relation coefficient) would confirm the existence of a random sec of vortices 
and would confirm the fluid dynamic model (Melick's analysis) representing 
the random sec of vortices. However, the data for the test conditions cor- 
responding to Figure 109 indicate chat the null hypothesis can be rejected 
with more than 99.95 percent probability, i.e.. Probability (derived cross 
correlation coefficient ^ measured cross correlation coefficient) ■> .9995. 
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However, the evaluation of these statistical data is inconclusive as to the 
determination of the existence or non-existence of a discrete set of random 
vortices. It is possible for the vortices to exist and the mathematical 
relationship between the vortex radius and cross correlation coefficient to 
be inaccurate, therefore yielding low probability of similarly between the 
measured and derived cross correlation coefficients. 

In summary, Melick’s procedure predicted peak distortion values to 
within 11.3 percent average error for fourteen data points having an average 
turbulence of 0.0155 and to within 20 percent average error for eight data 
points where engine stalls had been induced having an average turbulence of 
0.109. The procedure is inaccurate for test conditions with high levels of 
inlet turbulence. The derived and measured cross correlation data is incon- 
clusive as to the determination of the existence or non-existence of a dis- 
crete set of random vortices. If the vortices are present, then the theory 
for mathematical modeling of the vortices is inaccurate, and if the mathe- 
matical model of the vortices is correct, then the vortices are non-existent. 
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ENGINE STABILITY AUDIT EVALUATION 


P&WA has established a stability audit system for the F100 engine which 
has been successfully used throughout the F-15 development program. The capa- 
bility of the P&WA stability audit procedure to predict F100 engine stalls was 
evaluated for the ten flight conditions of Figure 110. Engine stalls were 
induced during the NASA F-15 Inlet/Engine Compatibility Flight Test Programs 
by positioning the inlet ramps to an off-design position, thereby artificially 
creating high inlet disrortiou levels. Ducing this test program, aircraft 
instrumentation was installed to accurately measure inlet distortion and 
pertinent engine parameters necessary to conduct stability audits. 

Stability Audit Procedure 

The F100 engine stability audit procedure was developed to predict fan 
and high compressor stability level for steady state operation from idle to 
maximum augmentation. In addition, it can also be used to predict fan and 
high compressor stability levels for normal augmentor transients within in- 
termediate power to maximum augmentation power setting. The engine audited 
in this study is an F100 (2 7/8) engine with Serial Number P690063. The F100 
(2 7/8) engine compression system components consist of an Improved Stability 
Fan (ISF) with a -25° to -5° inlet guide vane schedule and an Advanced Aero- 
dynamic Compressor (AAC) with a Scheme VII inlet guide vane schedule. 

An engine matching procedure was established which matched FI 00 (2 7/8) 
data obtained during sea level tests, flight tests and simulated ground test- 
ing at AEDC. Sufficient instrumentation was available during the NASA F-15 
Inlet/Engine Compatibility Test Program to obtain complete engine match points 
for the stability audits by matching the following measurands with the engine 
matching procedure. The procedure then provided predictions for any other 
required parameters. 

P C 2 ^ Fan inlet average total pressure 

T^ Fan inlet average total temperature 

N^ ''■> Low rotor speed 

CIW Fan inlet guide vane angle 

Pt 6 /P Cl ^ Engine pressure ratio 
m 

RCW ^ High Pressure Compressor (HPC) inlet guide vane angle 

The surge line destabilizing influences accounted for in the stability 
audit include Reynolds number, engine-to-engine variations and inlet distor- 
tion as shown in Figure 111. Operating line destabilizing influences accounted 
for include distortion rematch and augmentor combustion anomalies. 

The surge line destablizing influences applied to the highest available 
surge line (Point A) are shown in Figure 111. These influences include Rey- 
nolds number (Point B) , engine-to-engine variations (Point C) and inlet dis- 
tortion (Point D) . The operating line destablizing influences applied to 
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the original natch Point 1 are also shown. These include distortion rematch 
(Point 2), augmentor blowout (Point 3) and augmentor reignition (Point 4). 
After accounting for all destabilizing influences, the stall margin remaining 
is determined by subtracting the final operating pressure ratio from the final 
surge pressure ratio and normalizing this with the highest available surge 
pressure ratio, e.g. (PR^ - PR^/PR^. 

Surge line losses due to Reynolds number effects are based on a correla- 
tion of surge line loss versus Reynolds number for various engines including 
the F100. These losses were applied to fan and HPC undistorted surge lines 
previously obtained from numerous rig and engine tests. The surge line loss 
attributed to engine- co-engine variation was audited at one percent for both 
the fan and high compressor. 

The surge line loss due to inlet distortion was based on time variant 
distortion measured with the 48 probe engine face rake. Fan and HPC sensi- 
tivities to inlet distortion have been previously determined from F100 rig 
and engine tests. Peak time variant distortion values were determined using 
an analog distortion calculator to "flag” maximum time variant distortion 
levels for non-stall test conditions. For test conditions with engine stall, 
values occurring 8-20 milliseconds prior to the engine stall were used for 
the audits. 

The destabilizing effect on fan and high compressor operating lines due 
to distortion rematch was 1% and 0% respectively as determined from data 
collected during previous engine testing. 

For sore of the test conditions where engine stalls were inducad by 
fully extending the third ramp actuator an augmentor anomaly (blowout and 
reignition) would occur. The effect of this anomaly on the fan operating 
line has been determined from measured engine parameters and factored into 
the stability audits of this study. 

Induced Stall Procedure 

The data for the induced engine stalls were obtained on a F100(2 7/8) 
test engine with Serial Number P690063 which was mounted in the left hand 
side of Aircraft Number 71-0281. Stalls were induced by manually extending 
the left inlet third ramp actuator shown in Figure 2 beyond the scheduled 
position. 

Extending the inlet third ramp actuator beyond the scheduled position 
results in a smaller inlet throat area which significantly reduces total 
pressure recovery and increases fan distortion factor, as shown in 
Figure 112. Note the reduction in total pressure recoverv“and increase in 
distortion which occurs relative to normal third ramp scheduling (13 * 11°). 
This total pressure reduction brings about an increase in the effective 
pressure altitude at the engine face. The "as tested" altitude and Mach 
number for each of the stability audit points as well as the engine face 
effective pressure altitude is shown in Figure 113, 
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Stability Audit Test Conditions 

The ten test conditions of Figure 110 were selected to provide represen- 
tative stability audits throughout the F-15 flight envelope. Engine stalls 
were induced for eight of the ten test conditions by fully extending the third 
ramp actuator. The Mach number range for these eight test conditions is 0.4 
to 1.2. The remaining two test conditions are for Mach number 1.6 and 2.0. 

For flight safety reasons, engine stalls were not induced at the high super- 
sonic Mach number flight conditions. The instrumentation pertinent to mea- 
suring time variant fan distortion and for detecting the stalling component 
(fan and high pressure compressor) and augmentor instabilities is shown in 
Figure 114. 

During Che process of inducing engine stalls on engine P690063, augmentor 
combustion anomalies occurred for three of the ten test conditions. For Data 
Points 3, 4 and 16, augmentor blowouts occurred and were followed by stalls 
Induced by augmentor reignition pressure pulses. The augmentor blowouts 
occurred because the augmentor was operating outside the design envelope. 

This occurred because the repositioned inlets ramps reduced the engine face 
total pressure to a level corresponding to an altitude outside the design 
envelope, as shown in Figure 113. The stability audit can account for the 
effect of augmentor blowout and reignition on the fan operating point by using 
the fan discharge high response probe data. 

Upon examining the time history traces of fan distortion, it was found 
that the fan would, for some test conditions, tolerate distortion levels 
higher than those that finally stalled the fan. As the total pressure at the 
engine face was reduced, due to third ramp actuator extension, the engine 
operating pressures were reduced and the fan would match to a lower operating 
fan pressure ratio, thereby allowing tolerance to higher distortion levels. 

As the engine control system matched the fan back up to the normal operating 
fan pressure ratio, a peak distortion level would be encountered of sufficient 
level to stall the fan, thereby explaining the engine's temporary tolerance 
to higher than stalling distortion levels. These fan operating point excur- 
sions are depicted on a fan map in Figure 115. 

Stability Audits 

Five stability audits representative of the ten audits conducted are 
presented in the set of Figures 116 through 120. The audits for the remain- 
ing five test conditions are presented in Appendix J. The first figure of 
each set presents data traces of fan discharge high response pressures, fan 
inlet high response pressure and inlet third ramp actuator position at the 
time of the stall point audited or for the time at which peak distortion 
occurred during a non-stall poinc. The second figure of each set indicates 
the peak distortion value selected just prior to stall or the peak value 
during the event for a no-stall point. The third figure of each set presents 
the fan inlet pattern audited for each flight point, and the fourth figure 
of each set presents the fan and HPC stability audit results for each test 
condition. 

For the fan and HPC stability audits, the corrected airflows presented 
were converted from the engine value to an equivalent fan and HPC "rig basis". 
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Stability audits are conducted on this basis since the original stall margin 
loss and distortion correlations were established from rig data. Stability 
audits for five of the ten test conditions are discussed below. 

Stability Audit (Mach 0.59, Altitude«7610m, WAT2=104.1^ > I.D.»2) - A fan 
stall during steady state augmentor operation with the third ramp actuator 
fully extended is illustrated in Figure 116a. An audit of the peak distortion 
just prior to stall (see Figures 116b and 116c) indicate a negative 6.2 percent 
fan stall margin and a positive 18.4 percent HPC stall margin remaining as 
shown in Figure llod. 

Stability Audit (Mach 0.52, Altitude=10960m, WATZ^IO? .1Z, I.D.»3 ) - A 
rumble induced augmentor blowout followed by an augmentor auto-ignition re- 
light, which induced a fan stall with the third ramp actuator full extended 
is illustrated in Figure 117a. An audit of the peak distortion prior to surge 
from Figure 117b indicates stable fan operation after the blowout but prior 
to augmentor reignition. Once reignition occurs, a negative 27.4 percent fan 
stall margin is shown in Figure 117d along with 14.7 percent positive HPC stall 
margin remaining. The augmentor blowout and reignition effects on the fan 
operating point were determined from the fan discharge high response pressure 
probe traces. 

Stability Audit (Mach 0.35, Altitude =107fen, WAT2=104.2”, I.P.=15 ) - A 

fan stall during steady state augmentor operation with the third ramp actuator 
fully extended is illustrated in Figure 118a. Fan stability audits of the 
highest non-stall distortion, as well as the peak stall inducing distortion, 
from Figure 113b are presented in Figure USd. The non-stall case indicates 
-H3.6 percent positive fan stall margin remaining while the stall case indicates 
a negative 6.4 percent fan stall margin. The HPC stability audit is shown in 
Figure USe and indicates a positive 15.0 percent HPC stall margin remaining. 

Stability Audit (Mach 1.2, Altitude=li220m. WAT2=98.3”, I.O.-34 ) -A fan 
stall during steady state augmentor operation with the third ramp actuator 
fully extended is illustrated in Figure 119a. An audit of the peak distortion 
prior to stall from Figure 119b indicates a negative 3.1 percent fan stall 
margin and a positive 16.2 percent HPC stall margin remaining as shown in 
Figure 119d. 

Stability Audits (Mach 1.6, Altitude=18470in. *'AT2=39.3, I.D.=c4 ) - High 
response pressures and inlet third ramp actuator positon at the time of peak 
distortion are illustrated in Figure 120a for this test condition. This was 
a supersonic non-stall event with the engine power set at maximum augmentation, 
and the inlet third ramp automatically scheduled to its normal position. The 
peak distortion factor shown in Figure 120c was audited and the results are 
shown in Figure 120d. The audits indicate a positive 22.4 percent fan stall 
margin remaining and a positive 17.0 percent HPC stall margin remaining. 

Screening Curve Verification 

The levels of distortion values for which a stall may occur for the F100 
(2 7/3) engine are shown in the distortion screening curve of Figure 121. It 
is a screening curve in the sense that the values from this curve are not to 
be used as distortion limits for engine stall, but when an inlet distortion 
level for a particular airflow exceeds the comparable value from the screening 
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curve, additional analysis should be conducted to determine if there is any 
remaining engine stall margin for the particular test condition. The fan 
distortion levels for the audited test conditions are plotted on the screening 
curve of Figure 121. For the eight test conditions where engine stalls were 
induced, the distortion levels lie above the line in the region of predicted 
stability problems and the distortion levels for the two no stall test condi- 
tions lie below the screening curve in the region where stable operation is 
predicted. These results lend credibility to both the stability audit system 
and the screening curve. 

Identification of Stall Inducing Distortion Peaks 

For each of the stall points audited, the peak time variant distortion 
value just prior to the stall was chosen as the level which induced the stall. 
For each test condition evaluated, the data record time at which the peak dis- 
tortion and the engine stall occurred are summarized in Figure 122, and the 
elapsed time from the beginning of the distortion data record to the engine 
stall is presented. The elapsed time from the audited peak distortion to the 
stall point, as indicated by stall overpressure, ranged from 8 to 20 milli- 
seconds for all events. Numerous studies have indicated that this is the 
time interval magnitude required for the peak distortion pressure wave to 
travel from the instrumentation plane to the engine stalled compression stage, 
for the stall to occur, and for the stall overpressure to return to the instru- 
mentation plane. 

In summary, evaluation of the stability audit system indicates that suf- 
ficient technology exists to successfully predict inlet/engine compatibility 
status during the early phase of the inlet/engine development programs. 


38 



CONCLUSIONS AND RECOMMENDATIONS 


Conclusions 

The following conclusions can be drawn from the results of this pro- 
gram: 

o Peak time variant distortion values from subscale inlet model wind 
tunnel tests are representative of those from full scale flight 
test aircraft. 

o The time variant pressure data are random stationary data, thereby 
allowing valid statistical analyses to be conducted. 

o The effect of the engine presence on total pressure recovery, peak 
time variant distortion and turbulence level is small but favorable. 

o The Reynolds number/scale evaluation indicates a general trend of 
increasing total pressure recovery, decreasing peak time variant fan 
distortion and decreasing turbulence level as a function of increas- 
ing Reynolds number/ scale. 

o The frequency content evaluation indicates that peak time variant 
fan distortion and turbulence level increase as a function of in- 
creasing filter cutoff frequency for all of the data evaluated 
in this study. 

o The capability 7 of the Pratt & Whitney Aircraft stability audit system 
to predict engine stalls has been verified for both stall and non- 
stall flight test conditions. 

o Predictions of peak distortion values using Melick's procedure 

are accurate to 11.3 percent average error for fourteen data points 
having nominal turbulence levels and are accurate to 20 percent 
average error (the maximum error approaches AO percent) for eight 
data point having high turbulence levels. 

Recommendations 

Based on the results of this program the following recommendation is 

made: 

o A statistical study should be made to evaluate the properties of the 
time variant distortion data. Specifically the distribution and 
variance of the peak time variant distortion values should be 
investigated as a function of analysis time to further aid in the 
interpretation of the data for this study. 
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FIGURE 5 
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FIGURE 9 

MCAIR FLIGHT TEST DATA MANAGEMENT SYSTEM 
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FIGURE 10 
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FIGURE 1 0 (Continued) 

MCAIR WIND TUNNEL TEST DATA PROCESSING SYSTEM 
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FIGURE 15 (Continued) 
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FIGURE 16 
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Fan Circumferential Distortion Descriptor ~ Kg Fan Radial Distortion Descriptor ~ b K r a2 
Fan Distortion Descriptor ~ Kg + b K ra ^ High Compressor Distortion Descriptor ~ K C2 


^otniNAE 


59 





DISTORTION 

DESCRIPTOR 



PERCENT CORRECTED ENGINE AIRFLOW 


SUHHUIt 


FIGURE 17 

FAN DISTORTION SCREENING CURVE 


RAOIAL 

DISTORTION 

WEIGHTING 

FACTOR 

b 



FIGURE 18 

FI 00(2 7 8) RADIAL DISTORTION WEIGHTING FACTOR 


60 






FIGURE 19 

TEST CONDITIONS FOR EVALUATING 
QUALITY OF TIME VARIANT DATA 


s. 






















































awt wo n 


FIGURE 20 
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Random vs Periodic Data 
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FIGURE 24 

INSTANTANEOUS. MEAN AND MEAN SQUARE PRESSURE TIME HISTORY 

FLIGHT TEST 
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FREQUENCY - HERTZ 


Bl AUTOCORRELATION FUNCTION C) PROBABILITY DENSITY FUNCTION 


NORMALIZED 

AUTOCORRELATION 

FUNCTION 



LAG TIME SEC STANDARD DEVIATION 


FIGURE 25 

TIME VARIANT DATA QUALITY 
FULL SCALE WITH ENGINE 

Mach 1.8 a = — 2 (3 = 0 p=2.9 A 3 =18.6 

Bypass = 0 WAT2 = 79.8% CIW = —25.0 I D Number = 52 
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VIEW OF LEFT 
INLET LOOKING 
DOWNSTREAM 
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TIME -SEC 
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FIGURE 26 

INSTANTANEOUS. MEAN ANO MEAN SQUARE PRESSURE TIME HISTORY 

FULL SCALE WITH ENGINE 

Mach 1 8 a - —2 0 = 0 p = — 2.9 .A*=1R6 

Bypass = 0 WAT? -r 79 . 5 % ~; vv = —25.0 I D Number = 52 
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FIGURE 27 

TIME VARIANT DATA QUALITY 
FULL SCALE COLD PIPE 

Mach 1.8 a = — 2, 0=0 p = — 3.0 _W= 18 7 

Bypass = 0 WAT2 = 82.2% CiW = —25.0 I D Number = 49 
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VIEW Of LEFT 
INLET LOOKING 
DOWNSTREAM 


FIGURE 28 

INSTANTANEOUS. MEAN AND MEAN SQUARE PRESSU RE TIME HISTORY 
FULL SCALE COLD PIPE MODEL 

Mach 1.8 a = — 2 0=0 p = — 3.0 A 3 =18.7 
Bypass = 0 WAT2 = 82.2% CIW = —25 0 I D. Number = 49 
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FIGURE 29 

TIME VARIANT DATA QUALITY 
1 /6th SCALE INLET MODEL 

Mach 1.6 « = — 4 0 = 0 p = — 2 A 3 =13.5 

Bypass = 0.0 WAT 2 = 87.3% CIW = —21 .8 I D. Number = 42 
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FIGURE 30 

INSTANTANEOUS. MEAN ANO MEAN SQUARE PRESSURE TIME HISTORY 

1/ 6th SCALE INLET MODEL 

Mach 1.6 a = —4 0=0 p = — 2.0 A 3 =13.6 

Bypass = 0 WAT2 = 87.3% CIW = — 21.8 1.0. Number = 42 
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I INTERPOLATION 
I PROCEDURE 


SOLID SYr.’BOL STEADY STATE 
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TIME VARIANT FAN 
DISTORTION 


sym| 

MOD c *_ 1 

I D. NUMBER 

o 

1 6th SCALE 


□ 

FSCP 

49. 50 

O 

FSE 

51.52 

A 

FLT 

53 



FULL SCAL:: FlcTER CUTOFF FREQUENCY H* 

CP7MI2^}« 


FIGURE 31 

EFFECT OF FILTER CUTOFF FREQUENCY ON 
PEAK FAN DISTORTION PARAMETERS 

n---2 ft = 0 p = -3 A 3 =18 7 
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FIGURE 32 

DISTORTION DESCRIPTORS vs PERCENT CORRECTED ENGINE AIRFLOW 

Mach 8 a = — 2 (i - 0 p - — 3 A3 = 18 7 Bypass 0 
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FIGURE 33 

INTERPOLATION TECHNIQUE FOR FAN DISTORTION 
DESCRIPTORS vs PERCENT CORRECTED ENGINE AIRFLOW 

Mach 18 a - — 2 0-0 p - — 3 A 3 - 187 Bypass - C 
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Solid symbol - steady state 
Open symbol - corresponds to 
peak time var.ant tan distortion 

Flags indicate interpolated or 
extrapolated data relative to airflow 




FULL SCALE FILTER CUTOFF FREQUENCY Ha 
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FIGURE 34 

INTERPOLATION TECHNIQUE FOR FAN DISTORTION 
EFFECT OF FILTER CUTOFF FREQUENCY ON INSTANTANEOUS DISTORTION PARAMETERS 

Mach 1 8 a - -2 0 - 0 p= — 3 ^ 3 -1o7 EKpass 0 W4T2 - 78 9% 
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FIGURE 35 

TEST CONDITIONS TO EVALUATE THE EFFECT ON DISTORTION PARAMETERS 
OF "WITH" AND "WITHOUT" AN ENGINE 
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FIGURE 36 

COMPARISONS OF STEADY STATE AND 
INSTANTANEOUS TOTAL PRESSURE CONTOURS 
Mach 1.8 a = — 2 0=0 p = —3.0 

A 3 = 1 8.7 Bypass = 0 CIW = —25.0 Alt = 1 7.700 
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FIGURE 37 

COMPARISONS OF STEADY STATE AND 
INSTANTANEOUS TOTAL PRESSURE CONTOURS 
Mach 1 88 a = 4 0 =0 =2.5 

A 3 = 18.7 Bypass = 0 CIW = —25.0 Alt = 1 7.700 
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FIGURE 38 

COMPARISONS OF STEADY STATE ANO 
INSTANTANEOUS TOTAL PRESSURE CONTOURS 
Mach 2.2 a = — 2 S - 0 p = -4.0 
A 3 - 25.0/24.8 Bypass = 0.0278 0.0774 ClVV = -25 0 Alt - 19,800 - 21,300 
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FIGURE 39 

COMPARISONS OF STEADY STATE AND 
!NST ITANEOUS TOTAL PRESSURE CONTOURS 

•>2.2 o=4 0=0 p = 0.0/10 

A 3 = 2S.0 Bypass = 0.0774 CIW = —25.0 Alt * 1 9.800 - 21 ,300 
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FIGURE 40 

COMPARISONS OF STEADY STATE AND 
INSTANTANEOUS TOTAL PRESSURE CONTOURS 
Mach 2.2 a = 11 0 = 0 p = 6.8 

A 3 = 25.0/24.8 Bypass = 0.0677 - 0.0774 Cl VV = -25.0 Alt = 19,800 - 21,300 


81 




CN6M* 

Mfswr 

CfffCT 


MTHOUT BMGim iFSCPi 





47 A* * nr 3 
47.54 x HT 4 
47.54 x 10~® 
47A4X10 - * 
4734 x 10 -7 


MTH ENGINE (FSEI 


\ 


10. NUMBER • 52 
% WAT2 - 75.8 
Tu^7q » 0.0144 


47.54 * 10~* ' 



IfP 10 1 10 2 10 3 10 4 10® 10° 

FREQUENCY - Ht 


10 1 10 2 10 3 10* 
FREQUENCY ■ H* 


10 ® 


FIGURE 41 

COMPARISON OF POWER SPECTRAL DENSITY PLOTS 

Mach 1.8 a = — 2 0=0 p = — 3.0 A =ia7 CIW= — 25.0 
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FIGURE 42 

COMPARISON OF POWER SPECTRAL DENSITY PLOTS 

Mach 1.8 a = 4 0 = 0 p = 2.5 

A 3 = 1 8.7 CIW = —25.0 
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FIGURE 43 

COMPARISON OF POWER SPECTRAL DEN SI TY PLOTS 

Macn 2.2 a = — 2 00 p -- —4.0 A 3 = 25.0 CIW =' — 25 0 
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FIGURE 44 

COMPARISON OF POWER SPECTRAL DENSITY PLOTS 

Mach 2.2 a = 4 13 -0 p = 0 A 3 = 25.0 CIW = —25.0 
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FIGURE 45 

COMPARISON OF POWER SPECTRAL DENSITY PLOTS 

Mach 2 2 a = 11-12 0 = 0 p = 6 8 A 3 = 25 0 CIW = — 25 0 
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FIGURE 46 

TABULATED DATA COMPARISON OF AVERAGE HUB AND TIP STATIC PRESSURES FOR FULL 
SCALE INLET MODEL WITH AND WITHOUT ENGINE PRESENCE 



FIGURE 47 

COMPARISON OF AVERAGE HUB AND TIP STATIC PRESSURES FOR FULL SCALE INLET 

MODEL WITH AND WITHOUT ENGINE 
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FIGURE 48 

WITH AND WITHOUT ENGINE COMPARISONS OF DERIVED TURBULENCE 
LEVELS vs FILTER CUTOFF FREQUENCY 
Mach 18 « = — 2 0 = 0 p = —3 0 A 3 = *. 8 7 C!W = —25 0 
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FIGURE 49 

WITH AND WITHOUT ENGINE COMPARISONS OF DERIVED TURBULENCE 
LEVELS vs FILTER CUTOFF FREQUENCY 

Mach 1 8 « = 4 /i=0 = 2.5 ^ 3 =18 7 CIW = — 25.0 
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INLET LOOKING 
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Note: These Jerived turbulence levels are 

computed bv integrating the respective 
PSD plots. 
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FIGURE 50 
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WITH AND WITHOUT ENGINE COMPARISONS OF DERIVED TURBULENCE 


LEVELS vs FILTER CUTOFF FREQUENCY 


Mach 2.2 <» = — 2 0 0 p^-40 _* 3 = 25 0 CIW = -25 0 
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FIGURE 51 

WITH AND WITHOUT ENGINE COMPARISONS OF DERIVED TUR3ULENC r . 
LEVELS vs FILTER CUTOFF FREQUENCY 
Mach 2 2 a = 4 J - 0 p - 0 .A 3 = 25 0 CIVV - —25 0 
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FIGURE S3 

COMPARISON TOTAL PRESSURE RECOVERY 
WITH AND WITHOUT ENGINE 
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AVERAGE TURBULENCE WITHOUT ENGINE Tu 


PEAK FAN DISTORTION WITHOUT ENGINE K, 2 


*0 * open bypass. C * closed bypass. P - partially open bypass snasa 

FIGURE 54 

COMPARISON OF PEAK FAN DISTORTION AND TURBULENCE 
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*0 = Open Bypass. C = Closed Bvpass 

"These data exhibited harmonic charactenst.es and are not included in this analysis 
** 'These data have different bypass areas and are not included in this analysis 


FIGURE 55 

TEST CONDITIONS FOR EVALUATING REYNOLDS NUMBER SCALE EFFECTS 
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FIGURE 56 


REYNOLDS NUMBER AND MACH NUMBER RANGE OF THE FLIGHT TEST 
FULL SCALE AND 1/6th SCALE DATA 
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EFFECTS 


Flags indicate interpolated or 
extrapolated data relative to airflow 
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FIGURE 57 

EFFECT OF REYNOLDS NUMBER/SCALE ON TOTAL PRESSURE 
RECOVERY AND TURBULENCE 

Wacn0 6 = — 10 /3-10 p - —3 A 3 -10b Bypass = 0 0 
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FIGURE 58 

EFFECT OF REYNOLDS NUMBER/SCALE ON FAN 
DISTORTION AND SPATIAL DISTORTION 

Mach 0.6 a - — 10 /? - 10 p - —3 A 3 =106 Bypass 00 
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extrapolated data relative to airflow 
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FIGURE 59 

EFFECT OF REYNOLDS NUMBER/SCALE ON TOTAL PRESSURE 
RECOVERY AND TURBULENCE 

Mach 0 6 a = 4 p = 0 p-1 A3 = 10.6 Bypass = 0 0 
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FIGURE 60 

EFFECT OF REYNOLDS NUMBER SCALE ON FAN 
DISTORTION AND SPATIAL DISTORTION 

o=4 ;i - 0 .i = 7 A 3 = 10 6 Svpass -00 
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FIGURE 61 

EFFECT OF REYNOLOS NUMBER/SCALE ON TOTAL PRESSURE 
RECOVERY ANO TURBULENCE 

Mach 0 9 a - — 1C 0 = 10 p = 3 - 10.6 Bypass = 00 
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REYNOLDS 
NUMBER SCALE 
EFFECTS 


Solid symbol • steady state 
Open symbol - corresponds to 
peak time variant Lm distortion 

Flays indicate interpolated or 
extrapolated data relative to airflow 
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FIGURE 62 

EFFECT OF REYNOLDS NUMBER SCALE ON FAN 
DISTORTION AND SPATIAL DISTORTION 
Mach 0 9 « 1 —10 0 - 10 p - 3 ^3 = 10 6 Bypass = 00 
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extrapolated data relative to airflow 
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FIGURE 63 

EFFECT OF REYNOLDS NUMBER/SCALE ON TOTAL PRESSURE 
RECOVERY AND TURBULENCE 

Mach 0 9 a - — 4 0 = 0 P = — 1 A 3 = 8 2 


Bypass = 00 
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peak time variant fan distortion 

Flags indicate interpolated or 
extrapolated data relative to airflow 
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FIGURE 64 

EFFECT OF REYNOLDS NUMBER SCALE ON FAN 
DISTORTION AND SPATIAL DISTORTION 

Mach 0 9 u = -4 = 0 p = — 1 .13 = 82 Bypass = 0 0 
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REYNOLDS NUMBER x 10~ 6 BASED ON SCALE 


FIGURE 65 

EFFECT OF REYNOLDS NUMBER/SCALE ON TOTAL PRESSURE 
RECOVERY AND TURBULENCE 
Mach 0.9 a = 4 /3 = 0 p-1 A 3 = 10.4 Bypass = 00 
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Flags indicate interpolated or 
extrapolated data relative to airflow 
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FIGURE 66 

EFFECT OF REYNOLDS NUMBER/SCALE ON FAN 
DISTORTION AND SPATIAL DISTORTION 

Mach 0.9 a = 4 0=0 p -1 A3 - 104 Bypass = 00 
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extrapolated data relative to airflow 
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FIGURE 67 

EFFECT OF REYNOLDS NUMBER/SCALE ON TOTAL PRESSURE 
RECOVERY AND TURBULENCE 

Mach 1 2 a =10 0-0 p = 7 A 3 = 10 6 Bypass = 00 
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Open symbol - corresponds to 
peak time variant fan distortion 

Flags indicate interpolated or 
extrapolated data relative to airflow 
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FIGURE 68 

EFFECT OF REYNOLDS NUMBER/SCALE ON TOTAL PRESSURE 
RECOVERY ANO TURBULENCE 

Mach 1.2 or- 10 0=0 p= 7 ^3 = 10 6 Bypass = 00 


108 
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EFFECTS 


Flags indicate interpolated or 
extrapolated data relative to airflow 
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FIGURE 69 

EFFECT OF REYNOLDS NUMBER/SCALE ON TOTAL PRESSURE 
RECOVERY ANO TURBULENCE 

Mach 1 6 q = — 4 /3 = 0 p - —2 A3=13.5 Bypass = 00 
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Flags indicate interpolated or 
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FIGURE 70 

EFFECT OF REYNOLDS NUMBER /SCALE ON FAN 
DISTORTION AND SPATIAL DISTORTION 

Mach 1.6 a - — 4 j3=0 p=— 2 3 = 135 Bypass = 0.0 
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FIGURE 71 

EFFECT OF REYNOLOS NUMBER/SCALE ON TOTAL PRESSURE 
RECOVERY AND TURBULENCE 

Mach 1.8 a = -2 0 = 0 p = -3 .13=17.4 Bypass 0.9 
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FIGURE 72 

EFFECT OF REYNOLDS NUMBER/SCALE uN FAN 
DISTORTION AND SPATIAL DISTORTION 

Mach 1.8 a - —2 0 = 0 p = —3 A 3 = 1 7 4 Bypass = 0 0 
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FIGURE 73 

EFFECT OF REYNOLDS NUMBER /SCALE ON TOTAL PRESSURE 
RECOVERY AND TURBULENCE 

Mach 1.8 a = — 2 (3 = 0 p = — 3 .A 3 =18.7 Bypass = 0 0 
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Open symbol - corresponds to 
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FIGURE 74 

EFFECT OF REYNOLDS NUMBER/SCALE ON FAN 
DISTORTION AND SPATIAL DISTORTION 

Mach 1.8 or = — 7 0 ' 0 p - — 3 .A 3 - 18 7 Bypass = 0.0 
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Flags indicate interpolated or 
extrapolated data relative to airflow 
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FIGURE 75 

EFFECT OF REYNOLDS NUMBER/SCALE ON TOTAL PRESSURE 
RECOVERY AND TURBULENCE 

Mach 2.2 a = 0 3=0 p =— 2 ^ 3 = 22 5 Bypass = 0 0 
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FIGURE 76 

EFFECT OF REYNOLDS NUMBER /SCALE ON FAN 
DISTORTION ANO SPATIAL DISTORTION 

Mach 2 2 o = 0 P = 0 p=— 2 A 3 = 22.5 Bypass = 0.0 
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FIGURE 77 

EFFECT OF REYNOLDS NUMBER SCALE ON TOTAL PRESSURE 
RECOVERY AND TURBULENCE 

Mach 2 2 a - —2 0-0 p = — 4 .A 3 = 25 0 Bypass 0 071 
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FIGURE 78 

EFFECT OF REYNOLDS NUMBER/SCALE ON FAN 
DISTORTION AND SPATIAL DISTORTION 

Mach 2.2 a = —2 0 = 0 p = — 4 A 3 = 25 0 Bypass = 0 071 
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FIGURE 79 

EFFECT OF REYNOLDS NUMBER/SCALE ON TOTAL PRESSURE 
RECOVERY AND TURBULENCE 

Mach 2 2 cr=11 /3 = 0 p= 6 8 A 3 = 25 0 Bypass = 0 077 
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FIGURE 80 

EFFECT OF REYNOLDS NUMBER /SCALE ON FAN 
DISTORTION AND SPATIAL DISTORTION 

Mach 2 2 = 1 0 = 0 *> = 68 A 3 = 25 0 Bypass = 0 077 
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FIGURE 81 

EFFECT OF REYNOLDS NUMBER /SCALE ON TOT4L PRESSURE 
RECOVERY AND TURBULENCE 

Mach 2.5 er = 0 0 = 0 p = — 4 A 3 = 26.0 Bypass = 0 077 


121 






*«e blank 


*«T FUME, 


FREQUENCY 

CONTENT 

EFFECT 


Data Point 

Identification Number 




1/6* 

Scale 

Full 

Scale 

Cold 

Pipe 

Full 

Scale 

With 

Engine 

Flight 

Test 

*«o 

a 

(dig) 

4 

(deg) 

p 

(dag) 

Q 


Percent 

WAT2 

6 

- 

- 

7 

Eg 

-10 

10 

-3 

10.6 

C 

101.2 

5 

- 

- 


II 







17 

— 

— 

19 


-10 

10 

-3 

10.6 

C 

107.1 

18 

- 

- 


III 







42 

- 

— 

44 

1.6 


H 

-2 

13.5 

c 

89.1 

43 

- 

- 

- 



H 





- 

42,49.50 

51.52 

53 

1.8 

-2 

0 

-3 

18.7 

c 

78.9 

45 

— 

— 

47 

1.8 

-2 

n 

-3 

17.4 

c 

■2B 

46 

- 

- 




i 




1 

mrm 

62.63 


- 

2.2 

-2 

H 

-4 

25.0 

D 

KIM 

66 

68 

- 

70 

2.2 

m 

Hi 

-2 

22.5 

c 

73.0 

67 

69** 

- 



■i 

H 






79 

81 

— 

- 

2.5 

■ 

H 

mm 

26.0 

mm 

63.1 

80 

82 

- 



Hi 

H 

H 


H 



*0 * Open bypass. C * Closed bypass. P * Partially open bypass 

• 'These data exhibited harmonic characteristics and are not included in this analysis oan-aaas-rr 

'"These data ha«« different bypass door nosmons and are not included m this analysis 


FIGURE 83 

TEST CONDITIONS FOR 
EVALUATING FREQUENCY CONTENT EFFECT 
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FIGURE 84 

INTEGRATION OF THE POWER SPECTRAL DENSITY PLOT TO OBTAIN 
TURBULENCE AS A FUNCTION OF FILTER CUTOFF FREQUENCY 
Mach 2.2 a = — 2 0 = 0 p = —4 

A 3 = 25.0 WAT2 = 61 .7% Bypass = 0 I D Number = 62 
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TIME - SEC 


FIGURE 85 

COMPARISON OF TIME VARIANT FAN DISTORTION TIME HISTORIES 
FOR VARIOUS FILTER CUTOFF FREQUENCIES 

Mach 1.8 a = -2 0 = 0.0 p = -3 

A 3 = 18.7 WAT2 = 82.2% Bypass = 0.0 I.D. Number = 49 
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FIGURE 86 

EFFECT OF FILTER CUTOFF FREQUENCY ON FAN DISTORTION. 
TUf.AULENCE AND SPATIAL DISTORTION 

Mach 06 a - — 10 /3 = 10 p - — 3 A 3 =10.6 Bypass =0.0 
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FIGURE 86 (Continued) 

EFFECT OF FILTER CUTOFF FREQUENCY ON FAN DISTORTION, 
TURBULENCE AND SPATIAL DISTORTION 

Mach 0.6 a = — 10 - 10 p = —3 A 3 = 10.6 Bypass = 0 0 





FREQUENCY 

CONTENT 

EFFECT 


Solid symbol • steady state 

Open symbol - corresponds to Percent 

peak time variant fan distortion WAT2 

Flags indicate interpolated or ■ 107.1 

extrapolated data relative to airflow 


SYM 

MODEL 

I D. NUMBER 

0 

1/6th SCALE 

17,18 

□ 

FSCP 


0 

FSE 


A 

FLT 

19 


FAN DISTORTION 

K,, 


AVERAGE 

TURBULENCE 

<AP t J 
l 2 rms 


A) FAN DISTORTION 



O»7|0133M 


FIGURE 87 

EFFECT OF FILTER CUTOFF FREQUENCY ON FAN DISTORTION. 
TURBULENCE AND SPATIAL DISTORTION 

Mach 0 9 a = — 10 0=10 p - — 3 A 3 = 10.6 Bypass = 0.0 
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. IGURE 87 (Continued) 

EFFECT OF FILTER CUTOFF FREQUENCY ON FAN DISTORTION. 
TURBULENCE AND SPATIAL DISTORTION 
Mach 0.9 a - — 10 /3 * 10 p - — 3 A 3 = 10.6 Bypass = 0 0 
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FIGURE 89 (Continued) 

EFFECT OF FILTER CUTOFF FREQUENCY ON FAN DISTORTION. 

TURBULENCE AND SPATIAL DISTORTION 
Mach 18 a = —2 0=0 p = — 3 A 3 = 1 7 4 Bypass = 0 0 
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FIGURE 90 

EFFECT OF FILTER CUTOFF FREQUENCY ON FAN DISTORTION. 
TURBULENCE AND SPATIAL DISTORTION 

Mach 1.8 a - — 2 f$ = 0 q - — 3 -i 3 = 18.7 Bypass = 0.0 
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FIGURE 91 

EFFECT OF FILTER CUTOFF FREQUENCY ON FAN DISTORTION. 
TURBULENCE AND SPATIAL DISTORTION 

Mach 2.2 o = — 2 0 = 0 p = — 4 = 25 Bypass = 0 0774 
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FIGURE 92 

EFFECT OF FILTER CUTOFF FREQUENCY ON FAN DISTORTION. 
TUR 8 ULENCE AND SPATIAL DISTORTION 

Mach 2.2 a - 0 /3 - 0 j - —2 A 3 = 22.5 Bypass = 0.0 
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FIGURE 92 (Continued) 

EFFECT OF FILTER CUTOFF FREQUENCY ON FAN DISTORTION. 
TURBULENCE AND SPATIAL DISTORTION 

Mach 2.2 o=0 j}=0 p = —2 A3 = 22.5 Bypass = 0 0 
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FIGURE 93 

EFFECT OF FILTER CUTOFF FREQUENCY ON FAN DISTORTION. 
TURBULENCE AND SPATIAL DISTORTION 

Mach 2.5 a=0 0-0 p - —4 _A 3 =26 Bypass = 0 0774 
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Solid symbol - steady state 
Open symbol corresponds to 
peak time variant fan distortion 

Flags indicate interpolated or 
extrapolated data relative to airflow 
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FIGURE 93 (Continued) 

EFFECT OF FILTER CUTOFF FREQUENCY ON FAN DISTORTION. 
TURBULENCE AND SFATIAL DISTORTION 

Mach 2 5 <t = 0 (1=0 p = —4 .A3 r 26 Bypass = 0 0774 
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FIGURE 94 

EFFECT OF FILTER CUTOFF FREQUENCY AND REYNOLDS NUMBER 
ON FAN DISTORTION. TURBULENCE AND SPATIAL DISTORTION 

Mach 0.6 a = — 10 $ = 10 p = — 3 A 3 =10 6 Bypass = 0 0 
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FIGURE 94 (Continued) 

EFFECT OF FILTER CUTOFF FREQUENCY AND REYNOLDS NUMBER 
ON FAN DISTORTION. TURBULENCE AND SPATIAL DISTORTION 

Mach 0.6 a =—10 li = 10 p = — 3 A 3 =10.6 Bypass = 0 0 
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FIGURE 95 (Continued) 

EFFECT OF FILTER CUTOFF FREQUENCY AND REYNOLDS NUMBER 
ON FAN DISTORTION. TURBULENCE AND SPATIAL DISTORTION 

Mach 0 9 a =—10 0=10 p = — 3 A 3 = 10.6 Bypass = 0 
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FIGURE 96 

EFFECT OF FILTER CUTOFF FREQUENCY AND REYNOLDS NUMBER 
ON FAN DISTORTION, TURBULENCE AND SPATIAL DISTORTION 

Mach 1 6 a = — 4 /} = 0 p = — 2 A 3 =13.5 Bypass = 00 
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FIGURE 96 (Continued) 

EFFECT OF FILTER CUTOFF FREQUENCY AND REYNOLDS NUMBER 
ON FAN DISTORTION. TURBULENCE AND SPATIAL DISTORTION 

Mach 1.6 a-— 4 /? = 0 p = — 2 A3 =13.5 Bypass = 0.0 
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FIGURE 97 

EFFECT C r FILTER CUTOFF FREQUENCY AND REYNOLDS NUMBER 
ON FAN DISTORTION, TURBULENCE AND SPATIAL DISTORTION 

Mach 1 8 a = — 2 0 = 0 p - — 3 A 3 =17 4 Bypass = 0.0 
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FIGURE 97 (Continued) 

EFFECT OF FllTER CUTOFF FREQUENCY AND REYNOLDS NUMBER 
ON FAN DISTORTION. TURBULENCE ANO SPATIAL DISTORTION 

MacT. * 8 d = — 2 fi - 0 o~— 3 = 17.4 Bypass = 00 
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FIGURE 98 

EFFECT OF FILTER CUTOFF FREQUENCY AND REYNOLDS NUMBER 
ON FAN DISTORTION. TURBULENCE AND SPATIAL DISTORTION 
Mach 1.8 a - —2 j3 = 0 p - — 3 A 3 =18.7 Bypass = 00 
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FIGURE 99 

EFFECT OF FILTER CUTOFF FREQUENCY AND REYNOLDS NUMBER ON 
FAN DISTORTION. TURBULENCE ANO SPATIAL DISTORTION 

Mach 2 2 <> = -2 fi - 0 p = — 4 A 3 = 25 Bypass = 0 0774 
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FIGURE 100 

EFFECT OF FILTER CUTOFF FREQUENCY AND REYNOLDS NUMBER 
ON FAN DISTORTION. TURBULENCE AND SPATIAL DISTORTION 

Mach 2.2 o=0 /? = 0 p = 2 Jk 3 r 22.5 Bypass = 0.0 
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FIGURE 102 

TEST CONDITIONS FOR EVALUATING MEUCK S PROCEDURE 
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FIGURE 104 

COMPARISON OF MEASURED AND PREDICTED 
°FAK TIME VARIANT DISTORTION VALUES 
FOR INDUCED ENGINE STALLS 
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L4 

ENGINE FACE PRESSURE PROFILE 
LEFT FRONTAL VIEW 


Note: Inlet profile in terms of percent deviation from mean face pressure 

FIGURE 106 

PROBE SELECTION FOR CROSS CORRELATION COEFFICIENT CALCULATIONS 

Steady State Pressure Contour 

Mach 0.69 a = —8 4 fi - 10.6 I D. Number = 7 Flight 421 Run 10 
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‘Normalized cross correlation coefficient = 


[a xy It* * 0)1 


7 

9 

0 




“N.F. = Normalizing Factor 

FIGURE 109 

COMPARISON OF NORMALIZED VORTEX RADIUS AND CROSS CORRELATION 
COEFFICIENTS FROM THE ENGINE FACE RAKE PRESSURE PROBES 
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TEST CONDITIONS TO EVALUATE P8.WA STABILITY 
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Legend 
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B to C - Engine-to-engine variation 3 • Rematch due to augmentor blowout 
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FIGURE 112 

EFFECT OF THIRO RAMP ANGLE ON INLET PRESSURE RECOVERY AND FAN DISTORTION 
Mach 0.4 «* = *64 P = — 0.8 P = 69 WAT2 = 104 1 Bypass = 0.0 .O.Number>1 
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FIGURE 113 

COMPARISON OF F-15 FLIGHT ENVELOP? ANO TEST CONDITIONS AUOITEO 
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FIGURE 115 

FAN PRESSURE RATIO OPERATING POINT EXCURSIONS 
DURING THIRD RAMP ACTUATOR EXTENSION 
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FIGURE 116 

STABILITY AUDIT ANALYSIS PLOTS 

Mach 0.59 o = 13 9 0 = 0,9 #? = ? 0 

A 3 * 26.6 WAT2 * 102 7 % Bypass =■ 0 1 0 Number = 2 
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FIGURE 116 (Continued) 

STABILITY AUDIT ANALYSIS PLOTS 
Mach 0.5° a = 13.9 0 = 0.9 p = 7.0 

^3 = 26.6 WAT 102.7% Bypass = 0 1.0 Number = 2 
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FIGURE 116 (Continued) 

STABILITY AUOIT ANALYSIS PLOTS 
Mach 0 59 <* = 13 9 0=0 9 p = 7 

A 3 = 26 6 WAT2 = 102.7% Bypass = 0 I D Numb, r = 2 
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FIGURE 117 

STABILITY AUDIT ANALYSIS PLOTS 

Mach 0 52 a = 10.0 0 = 0.07 p = 7.0 
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FIGURE 117 (Continued) 

STABILITY AUDIT ANALYSIS PLOTS 

Mach 0.52 a = 10 0 0=07 p = 7 0 

27.6 WAT2 = 107 1% Bypass = 0 I D. Number = 3 
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FIGURE 117 (Continued) 

STABILITY AUDIT ANALYSIS PLOTS 

Mach 0.52 a = 10.0 0 = 0 07 p = 7 0 

A 3 = 27.6 WAT2 = 107.1% Bypass = 0 ID Number = 3 
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FIGURE 118 

STABILITY AUDIT ANALYSIS PLOTS 

Mach 0.85 a = 8.8 j3 ~ — 0.5 p = 7.0 

A3 = 27.6 WA72 - 104.2% Bypass = 0 I.D. Number =15 
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FIGURE 119 

STABILITY AUDIT ANALYSIS PLOTS 
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FIGURE 119 (Continued) 

STABILITY AUDIT ANALYSIS PLOTS 
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stability 

AUDITS 


Cl TOTAL PRESSURE CONTOUR AT PEAK FAN DISTORTION 

L-8 | TOP 



Lj£ 


FIGURE 110 (Continued) 

STABILITY AUDIT ANALYSIS PLOTS 

Mach 1.21 a =1.5 0=0 p = 60 

A 3 - 27 6 WAT2 = 98 3% Bypass = 0 1.0. Number = 34 
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STABILITY AUDIT ANALYSIS PLOTS 
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appendix a 


COMPARISONS OF DIGITIZED AND ANALOG DISTORTION DATA 

Presented herein are the comparisons of time history fan distortion 
values computed by the Analog Distortion Calculator and by the Digital 
Process of Volume I, Figure 10. A list of illustrations as provided below. 
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FIGURE A-2 

COMPARISONS OF DIGITIZED ANO ANALOG DISTORTION OATA 

M 0 = 0.59 a = 13.9 0 = 0.9 p- 7 0 

A 3 = 26 6 WAT2 = 102.7% Bypass = 0.0 
I D Number = 2 fc/o=170 
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FIGURE A-4 

COMPARISONS OF DIGITIZED AND ANALOG DISTORTION DATA 

M c = 0.69 o = 11.5 /? = 10 p = 7.0 
= 26 5 WAT2 = 104.2% Bypass = 0 0 
I D. Number = 4 fc. o = 170 
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FIGURE A-5 

COMPARISONS OF DIGITIZED AND ANALOG DISTORTION DATA 

M 0 = 0.69 a = —84 0 = 10.6 p= 0.6 

Jk 3 = 10.5 WAT2 = 101.2% Bypass = 0.0 
1.0. Number = 7 fc/o = 170 
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FIGURE A-6 

COMPARISONS OF DIGITIZED AND ANALOG DISTORTION DATA 

M 0 = 0.67 a = 4.3 0 = 0.7 P- 6.9 

A 3 = 11.1 WAT2 = 94 4% Bypass = 00 
I D. Number = 11 fc/o = 170 
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A) DIGITIZED 
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FIGURE A-7 

COMPARISONS OF DIGITIZED AND ANALOG DISTORTION DATA 

M 0 = 0.69 a = 3.4 $ = 0.7 p = 6 9 

A 3 = 11.1 WAT2 = 74.1% Bypass = 0.0 
I.D. Number =12 fc/o = 170 
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FIGURE A-8 

COMPARISONS OF DIGITIZED AND ANALOG DISTORTION DATA 

M 0 = 0.59 a= 4.6 ,8=1.2 p=7.0 

A 3 = 11 1 WAT2 = 107.9% Bypass = 0.0 
I D. Number = 1 3 fc/o = 1 70 
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FIGURE A*9 

COMPARISONS OF DIGITIZED AND ANALOG DISTORTION OATA 

M 0 = 0.85 a = 8.8 0=—O5 p= 70 

A 3 = 27.6 WAT2 = 1 04.2% Bypass = 0.0 
I D. Number = 15 fc/o=170 
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FIGURE A-10 

COMPARISONS OF DIGITIZED AND ANALOG DISTORTION DATA 

M 0 = 0 92 a = 5 6 0 = 0.6 p = 7.0 

^ 3 = 26.6 WAT2 = 104.5% Bypass = 00 
I D. Number =16 fc/o = 1 70 
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FIGURE A-1 1 

COMPARISONS OF DIGITIZED AND ANALOG DISTORTION DATA 

M 0 = 0.94 a = — 8.9 0 = 10.2 p=1.0 

A 3 = 10.5 WAT2 = 107.1% Bypass = 0.0 
I.D. Number = 19 fc/o = 170 
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FIGURE A-12 

COMPARISONS OF DIGITIZED AND ANALOG DISTORTION DATA 

M 0 = 0.90 a = —2.8 0 - —0.2 P = — 1 .2 

A 3 = 8.7 WAT2 = 97_5% Bypass * 0.0 
1.0. Number = 21 fc/o = 1 70 
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FIGURE A-13 

COMPARISONS OF DIGITIZED AND ANALOG DISTORTION DATA 

M 0 = 1.21 <* = 15 0 = 00 p = 60 

A 3 = 27 6 WAT2 = 98.3% Bypass = 0 0 
I.D. Number = 34 fc/o =170 
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FIGURE A-14 

COMPARISONS OF DIGITIZED AND ANALOG DISTORTION DATA 
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FIGURE A- 15 

COMPARISONS OF DIGITIZED AND ANALOG DISTORTION DATA 

M 0 =154 o=1.5 fi* 0.0 p=— 1.4 

A 3 = 27 0 WAT2 = 95 4% Bypass = 0.0 
1.0 Number = 41 fc/o = 170 
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FIGURE A- 16 

COMPARISONS OF DIGITIZED AND ANALOG DISTORTION DATA 
Mo* 1-75 a*— 2.6 0 = 0.4 p = -2 2 

A 3 = 1 6. 7 WAT2 = 80.7% Bypass = 00 
I D Number = 47 fc/o = 170 
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FIGURE A-17 

COMPARISONS OF OIGITIZEO AND ANALOG DISTORTION DATA . 

M 0 = 2.2 or = 0.3 /8 = 0.2 p = -2.2 

A 3 = 22.9 WAT2 = 73.0% Bypass = 0.0 
1.0. Number = 70 fc/o = 170 
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APPENDIX B 


DISTORTION DESCRIPTOR EQUATIONS AND SAMPLE CALCULATION 

Presented herein are the details of the distortion descriptor equations, 
of the o.stortion descriptor calculation procedure, of a sample calculation 
and of the distortion descriptor reference radial profile and radial distor- 
tion weighting factor that are used in this study. These topics are separated 
into sections as follows. 


Topic Page 


Summary of Distortion Descriptor Equations (see Figure B-l) 211 
Calculation Procedure 212 
Sample Calculation 216 
Reference Radial Profile & "b” factor (see Figure B-9) 233 
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FIGURE B-1 

SUMMARY OF DISTORTION DESCRIPTOR EQUATIONS FOR THE 
NASA/F-15 DISTORTION METHODOLOGIES STUDY 
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CALCULATION PROCEDURE FOR K a? AND K C2 FOR INL ET 
TIME VARIANT PRESSURE DATA 


1. Select the recorded inlet pressure pattern for the operating condition 
under consideration. 

2. This pattern is obtained from total pressure data measured at equidistant 
circumferential stations, 9^, $2, 83... etc., and at the mean of J annu- 
lar rings, dj_, d2, d3»..etc. , counting from the hub to the tip. 

3. Consider the ring d^ and list the total pressure ratio P t 2/P to for each 
angle of 9. 

4. For each angle 9 obtain the sine and cosine of the angle and the product 
(P t 2/ p to^ sine 9 and (Pt2/^to) cosine 9. 

5. Sum the values of Pt2/ p to sine 6 and P t 2/P to cosine 9 for all angles of 
9 to give 

TT IT 

£ P t2/ p to sine 9 30(1 T. P t2 /P to cosine 9 

~TT 

6. From the sum computed in step (5) and the value (A6/tt) , obtain the product 

1 T TT 

"T ^ P t2/ p to sine 9A9 and ^ £ p t2/ p to cosine 9A6 

T ' IT 

These values are respectively equal to a ^/ p t0 and b^/P^ 

7. The Fourier coefficient parameter Aj/ p to is then obtained 

Vto- V<n/f co ) 2 + <y V 2 

8. The value of the Fourier coefficient is determined for ten values of N 
where N (i.e., 9, 29, 39, ...108) is substituted for 9 in steps 4-7 and 
the largest value of 

is obtained. 

ring average value of pressure, P -*/P 

to 
Nq 

-i« 9 £ e t ,/P to 

is the number of circumferential location. 

10. Repeat steps 3-9 for each ring. 




P 2 
to N 


9. Find the 


P 0 */P 
t2* to 


where N, 


9 
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11. Multiply the values of l/P to (Ah/n 2) maximum obtained in step 8 by the 
weighting function 1/D and sum for all rings to obtain 

J 

L l/P to (A^/N^maximum 1/D 

Uhere J is the number of diameters. 

12. Obtain the sum of the weighting factors 1/D 

Li/d 

1 

13. Obtain P_,/P , the arithmetic average value of P^,*/P,. , from step 9 

CO CZ* CO 


P t2 /P to 


) 


L p t2* 


/p 


to 


14. Obtain Q/P^, where Q/P £ is defined as a function of percent total engine 
air flov corrected to trie engine face. 


15. Obtain Kg from the results of steps 11, 12, 13 and 14. 


Kq 


i 1 H 

|y. 


max 


(!ii) f-a.) 

V P to/avs VW 


J 1 
T-i- 

n 


t2' avg 


16. Obtain as from Q/P^ (step 14) and Figure B-2. 

17. Compute Splitter Diameter 


SPLITTER 


= Vas (° DZ ‘ 1dZ > + ID" 


18. Computer Kg as in step 15, but using values only for rings hav- 

ing diameter smaller than or equal to splitter diameter. 


19. Compute 


(K) 

\ t2*/ 


t2 ' inst 


(h) ’fa) M 

V t2*/inst \ to / inst /\ to /inst 


P 


213 



20. Plot 


vs 9 for ring diameters smaller than or equal to 


splitter diameter. 

21. Obtain 9~ where 9 Is equal to the greatest circumferential extent where 
Pt2/ p t2* <1.0. If there are two regions of low Pt2^ p t2* b y 25 deg or 
less, they are to be treated as one low pressure region. The lower limit 
of 9“ is to be 90 deg. 

22. Compute 


23. Compute 


c 2 “ K 0splitter 

- ■ feL/ft), 


24. Select appropriate reference radial profile. 

P - P 

25. For each ring, compute t2* t2* base 


P t2* ‘ P t2* base 


inst (step 23) 


_ p t2* 

p t2 base 


26. Compute 1 for each ring. 
D 


27. Comoute K 


P t2* " P t2* base 


± I zz 

£ (~l\ 

l t2 / avg 1 l D J 


t2 base/ P t2j 


28. Select appropriate radial distortion weighting factor 

29. Compute K_ 


Kq + b K 
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RATIO OF AVERAGE PRESSURE 
HEAD TO INLET TOTAL PRESSURE. Q/P t2 




FIGURE B-2 

F10G-PW-1C0 TURBOFAN ENGINE NOMINAL aS AS A FUNCTION OF THE RATIO 
OF AVERAGE DYNAMIC PRESSURE HEAD TO INLET TOTAL PRESSURE 
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Sample Calculation of and Distortion Descriptors 


Given Values: 



Diameter 


14.64036 

19.91952 

24.03636 

27.54072 

30.64836 

33.48108 



.9622 

.9783 

1.0016 

1.0166 

1.0191 

1.0185 


"b" factor 




See following page for 
probe pressure values 


( 7 ) 315 


Inboard 


( 6 ) 270 ' 


( 5 ) 225 * 



<> 3 

Rings 

< ) 2 


Outboard 



( 4 ) 

Looking Aft 


90 ° ( 2 ) 



135 ° ( 3 ) 


This leg is zenerally labeled (8). The distortion value 
computed with (8) is identical to this calculation. 

6 











Calculation ol Pfj 


Dlaaatar'S^ 

0 

45 

90 

135 

180 

225 

220 

315 

14.64036 

10.799600 

10.848300 

10. 89 .o99 

10.907100 

10.808999 

10.461499 

10.586100 

10.506100 

19.91952 

10.423499 

10.484099 

10.729799 

10.96999 

10.979699 

10.379700 

10.44100 

10. 260699 

24.03636 

10.183800 

10.232400 

10.777900 

10.690599 

10.733899 

10.368899 

10.239799 

10.143800 

27.54072 

10.099500 

10.231199 

10.868299 

10.629399 

10.550099 

10.275999 

10.511700 

9.9966993 

30.64836 

10.104600 

10.179600 

10.79400 

10.535399 

10.484200 

10.382.00 

10.176900 

9.5592991 

33.48103 

10.063200 

10.274799 

10.676900 

10.375600 

10.519099 

10.286799 

10.333400 

9.9924994 


*» <*t 2 > 

L. -Ts 1 ’ 10 * 46148 
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38 - 45*; 

M-l 


King 

■a 


Ring 81 


10.7496 


2.6999 


10.4235 




10.1838 




10.0995 



, 5 

10.1046 



pi IH 

10.0632 

Hi 



10.3483 

.1767767 


2 

10.4341 



3 

10.2324 



4 

10.2312 



5 

10.1796 



6 

10.2748 



1 

10.8919 

0.000^ 

0.0 

2 

10.7298 



3 

10.7779 



4 

10.8683 



5 

10.7940 



6 

10.6769 



1 

10.9071 

-0.1*0/777 

-1.92813 

2 

10.9700 



3 

10.6906 



4 

10.6294 



5 

10.5354 



6 

10.3756 



1 

10.8090 

-0.2SQ00Q 

-2.70225 

2 

10.9797 



3 

10.7339 



4 

10.5501 



5 

10.4842 



6 

10.5191 



L 

10.4615 

-0.176777 

-1.34936 

2 

10.3797 



3 

10.3698 



U 

10.2760 



5 

10.3827 



6 

10.2068 



1 

10.5861 

mmmmm 

0.0 

7 

10.4410 



3 

10.2398 

I 


4 

10.5117 



1 J 

10.1769 




10.3334 

I ■ 


, i 

'0.5061 

0.1767777 

1.85724 

1 2 

10. 2607 



3 

10.1.38 




9.9967 



J 

9.9593 



6 

9. 9925 



m 


-.00486 


I 
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Calculation of (a^) 


Ring * ^ c 2 


tP *2>rl«« * l5o C0 ‘ (Kk£ ®> 

Ring #2 

Ring *3 

ling *4 

Ring #5 

Ring #6 

mg u> 

2.605875 

2.54595 



2.5158 

0 


1 .30886 

g.. 80865 

1.79953 

1.31636 

1 

o.a 

m 

0.0 

0.0 

0.0 

2 

-1.93925 

— 

-1.38986 

-1.87904 

-1.86242 

-1.33417 

3 


Q 

-2.63752 

-2.62105 

-2.62978 

4 


■ 

-1.31657 

-1.83543 

-1.81348 

5 

0.0 

0.0 

0.0 

0.0 

0.0 

6 

■ 

1.79320 

1.76719 

1.76058 

■ 

7 

-.21598 

-.25848 

-.23242 

-.23264 

m 

^l’sing * ?c 2 
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1 
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1 

10.5001 

10.4410 
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-2.04052 


Hj 

-2.02792 


-2.581)5 

0 

1 

10.5001 

10.2007 

lo.uia 

9.9907 

9.959) 

9.9925 

-0.1707777 

-1.85724 

-1.81)80 

■ 

-1.70719 


-1 . 7664* 

7 
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10.4313 

10.3797 

10.3393 

10.2730 

10.3327 

10.2333 

-.173777 

-1.3494 

-1.3349 

-1.3331 

-1.8164 

-1.3354 

-1.3135 

5 

1 

10.3331 

10.4310 

10.2393 

10.3117 

10.1739 

10.3334 

0.2500 

2. 3435 

2.3102 

2.5300 

2.3279 

2.5442 

2.5833 

5 

1 

2 

3 

4 

s 

« 

10.3031 

10.2307 

10.1433 

9.9937 

9.9593 

9.9923 

-.173777 

-1.3572 

-1.3139 

-1.7932 

-1.7572 

-1.7606 

-1.7663 

7 

191 

1S33B9I 

.0323 

.0715 

-.0621 

.0143 

-.0835 

-.0203 

' ‘Vua, * f '2 
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Calculation of (fc ) 


* ft 


2 



A * 

- 43*; »•* 


< Vm »8 * IK 

■la 0*44) 



lln * 

■a 

usa 

llH *1 

tin * 42 

Ho * 43 

tint 44 

tin * #3 

81 a * 44 

La * (8) 

1 

io . m * 

0.000 

0.3 






0 

I 

10.4233 



0.0 






3 

10.1*3* 




0.0 





4 

10.0393 





0.0 




j 

10.104* 






0.0 



* 

10.0432 







0.0 


i 

10.04*3 

0.000 

0.0 






1 

2 

10.4041 



0.0 






3 

10.2324 




0.0 





4 

10.2312 





0.0 




3 

10.179* 






0.0 



« 

10.2740 







0.0 


1 

10.091* 

0.000 

0.0 






2 

2 

10.7290 



3.0 






3 

10.7779 




0.0 





4 

10.0403 





0.0 




5 

10.7940 






0.0 



4 

10.4749 



- 




3.0 



10.9071 

0.000 

0.0 






3 


10.9700 



0.0 







10.4904 




0.0 






10.4294 





0.0 




3 

10.5334 






0.0 



4 

10.3734 







0.0 



10.8090 

0.000 

0.0 






4 


10.9797 



0.0 







10.7339 




0.0 






10.3301 





0.0 




3 

10.4042 









! * 

10.3191 







0.0 



10.4413 

0.000 

0.0 






3 


10.3797 



0.0 







10.3494 




0.0 






10.2740 





0.0 




3 

10.3027 






0.0 



* 

10.2840 





. 


0.0 



10.3041 

0.000 

0.0 






6 


10.4410 



0.0 







10.2394 




0.0 






10.5117 





0.0 




5 

10.1/49 






0.0 



4 

10.3334 







0.0 



10.30*1 

1.000 

0.0 






7 


10.2607 



0.0 







10.1434 




0.0 






9.99*7 





0.0 




3 

9.9593 






0.0 



6 

9.9923 









0.0 


[ £ &>''2 * to 

< iw ,, ]* l «, * 

m 

m 

m 

0.0 

JEM 

mm 

ikflEUS 
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Stair of 4,/M 2 


«*« 

1 

2 

3 

4 

5 

6 

T 


r 

05 p t 2 oaOKM) 

-.OOAM 

•.2459* 

-.25*48 

-.23242 

-.23264 

-.1*3*2 

□ 

( Vria« * *»2 

t 

TO '*2 sia(KM) 

.2X373 

.21*05 

.20*9 

.19310 

.22022 

.15147 

D 

^ria* ' * fc 2 

t 

ns n 2 Vco. 2 («a»}««i* 2 c*a») 

.21571 

* 

.32739 

.33109 

.30217 

.32034 

.23*19 

fl 

*1*'«2 



1.2157* 

.32739 

.33109 

.30217 

.32034 

ec 

0 

*1 • ?c 2 /* 2 



■■■■■ 



■■■■■■■■ 


mmmJ 


i 

T5o't2 —<■“•) 

.012* 

.05*2 

-.0250 


-.0955 

-.1070 

2 

**2^ria* ' '*2 

t 

Yio *tj alaOBI*) 

-.025* 

-.091* 

-.05*1 

-.029* 

.017 

.04*4 

2 

^rlm* ' Pt 2 

t 


.04157 

.1087 

.06325 

.1850 

.0970 

.1174 

2 

4j * ?t 2 


■Ml 

.01039 

.02717 

.015*1 

.04625 

.02425 

.02936 

2 

4j * ft,/* 2 

i 

So r t z co*(«4») 

.0004 

-.0321 

-.0166 

.0072 

.0429 

-.044 

3 

‘Vrla, * ? -2 

- 

Ho P c 2 .ia(*4») 

.0626 

.0715 

-.0621 

.0146 

-.0*85 

-.0203 

3 

< b 3>rtM* ; *2 

T 

P t j ^coo^OKAi )4*la Z (JK£9) 

.0626 

.0784 

.06428 

.01628 

.09*35 

.04846 

3 

*3 * ? *2 



.00696 

.00871 

.00714 

.001*1 

.01093 

.0053* 

3 

*3 • 

: 

Ho F t2 coo(IK49) 

. .0*76 

.1198 

.1249 

.2240 

.1257 

.165* 


( Vrfn« ' '*2 

- 

Hq ? t2 siu(Mit) 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

B 

(fc 4 ) rla * * c 2 

z 

Hg- r C2 >|coo 2 (K4t)-*sla Z (IKAt) 

.0*7* 

.119* 

.1249 

.2240 

.1257 

.165* 

1 

*4. * '*2 



.00549 



.00749 

.00781 

.014 

.007*6 

.0104 

B 



* H*zlca A^M 2 Values 


Calculation of Kg 



1 

10.46148 


.21578 

14.64036 


.32739 .33109 .30217 

19.91952 24.03636 27.54072 


+ 


.32034 

30.64836 


.23819 

33.48108 


.061522 (.2589156) 


Kg - .44099 


TXCB* 

5‘ioo*° ,T, '" ry 
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DK POOR QUAUia 

Calculation of K r> 


'V 8 
IllnjS. 

0* 

45* 

90* 

135* 

180* 

225* 

270* 

315* 

King 

Average 


10.7996 

10.8483 

10.8919 

10.9071 

10.8090 

10.4615 

10.5861 

10.5061 

10.7261875 


10.42350 

10.4841 

10.7298 

10.9700 

10.9797 

10.3797 

10.4410 

10.26C7 

10.5835625 


10.1838 

10.2324 

10.7779 

10.6906 

10.7339 

10.3689 

10.2398 

10.1438 

10.4213875 


10.09950 

10.2312 

10.8683 

10.6294 

10.5501 

10.2760 

10.5117 

9.9967 

10.3953625 

5 1 

10.10460 

10.17960 

10.7940 

10.5354 

10.4842 ! 

10.3827 

10.1769 

9.9593 

10.3270875 

6 

10.06320 

10.2748 

10.6769 

10.3756 

10.5191 

10.2868 

10.3334 

9.9925 

10.3152875 


Avg P t j - 10.46144 



© 

© 

<8> 

©/© 

© 

Base Radial Profile 

© 

©-© 

© 

©/© 



— 


p c 2 */ p t0 


1 

AP 


A p t 2 

Ring 

Pa- * 
r to 

p i 

l 2 


"V>B.s. /F ‘2 

I 



P~ 





p t 2 ' p t 0 


1 

p t 2 


c 2 Ring 

M| 

10.7261875 

10.46148 

1.0253031 

.9622 

.0631031 

.0655821 

.•"Ml 

10.5835625 



1.0116697 

.9783 

.0333697 

.0341099 

Mi 

10.4213875 



.9961676 

1.0016 

.00543239 

.00542371 

Mi 

10.3953624 



.9936799 

1.0166 

.02292009 

.02254583 

MS 

10.3270875 



.9871535 

1.0191 

.03194641 

.03134767 

■ 

10.3152875 



.9860256 

1.0185 

.03247436 
— 

.03188450 




.0655821 .0341099 .00542371 .02254583 

" 14.64036 19.91952 24.03636 27.54072 


.03134767 

30.64836 


+ 


.03188450 

33.48108 


.061522 (.2589156) 


- .57827 






































Calculation of IL, 
*2 

Ka 2 * Kq + b * R r az 

- .44099 + (12.50) (.57827) 

K a2 - 7.6694 


Calculation of 


D * 

Splitter 


yji s (OD 2 - ID 2 ) + ID 2 


a * .47 from Figure B-2 
s 


yj.^7 (33. 48108 2 • 14.64036 2 ) + 14.64036 2 


D c . - 25. 33* *5 

Splitter 


Kco - Ke 


9 * 218 from Figure B-4 

180- 


Splitter 


6 


Ke 


Splitter 


£ 

Ring-1 


(?) 1 • 5 &)L * 2 2 U 

= d * ±9*3 J 


(\ ■ K\ 

1 

\ 2 / 

X D 

\ N /max 

Ring Ring 


[ (<i/ ^2» i £ 5-7- 

*- -J Ring-1 Ring 


[<q/? t 2 >] £ 

*“ -I Ring-1 


°Ring 


Ke, 


1 / . 21578 .32739 .33109 ) 

10.57705 \ 14.64036 19.91952 24.03636/ 180 


.061522 


(—1 + —1 + 

(14.64036 19.91952 24.036361 


218 


180 


Kc, - .43143 x - .35622 


Kc ? - .35622 


*This value usually computes so as to use the inner three rings. 


228 



Detemination of 9 


Maximum value of 9 occurs in Ring 2, which can be determined by examin- 
ing the Ring 1 through Ring 6 Circumferential Pressure Contours that are 
shown in Figures B-3 through B-8. In addition, the Fourier Series curve fit 
of Ring 2 Circumferential Pressure contour is plotted on the graph and the 
equation is shown below. The fit has less than .62 error throughout the 
entire contour. 

P 2 (9) - 10.5836 - .24598 cos (45K) + .21605 sin (45K) 

+ .0582 cos (90K) - .0918 sin (90K) 

- .0321 cos (135K) + .0715 sin (135K) 

+ .1198 cos (180K) 

K * 0, 1, 2, 3, 4, 5, 6, 7 
9 * 218° on Ring 2, Figure B-4 
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INDIVIDUAL 

PROBE 

PRESSURE 

PSIA 



FIGURE 8-3 

RING 1 CIRCUMFERENTIAL PRESSURE CONTOUR 


INDIVIDUAL 

PROBE 

PRESSURE 

PSIA 



FIGURE B-4 

RING 2 CIRCUMFERENTIAL PRESSURE CONTOUR 
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'Prototype F100 (2 7 8) Engine 


FIGURE 8-9 

REFERENCE RADIAL PROFILES AND RADIAL WEIGHTING FACTORS FOR THE 
NASA /F 15 DISTORTION METHODOLOGIES STUDY* 











Svabol 


SYMBOLS, SUPERSCRIPTS, AND SUBSCRIPTS 

Quantity 



H th coefficient of the sine ten in a 
* Fourier expansion of the p t2^ P to ** 

8 curve for one ring. 






* N C ** coefficient of the cosine term in 
a Fourier expansion of the P ./P^ vs 

A £ • Vfc tO 

8 curve for one ring. 


d. 

i 


* ring at diameter 


0 • diaaeter of ring - any uni ts 

f/f * normalized frequencv ratio 

a 

ID * Inside Dia*ueter 


J 


« Number of diameters 


K 

N 

■J 

OD 

P t 

b 


q 


* distortion factor 

* order of Fouriei coefficient (number 
of waves per circumference) 

* number of circumferential locations 

* Outside Diameter 

2 

* total or stagnation pressure - Ib/ft 

* radial distrotion weight factor 

» average dynamic pressure heat at com- 
pressor face - lb/ft* 


AQ * angular increment between two adjacent 

circumferential points. (This quanti- 
ty must be the same for all points.) 

- radians 
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AP 


c 


6 


T 

Superscripts 


(— ) 


( >* 

Subscripts 


avg 

c 

i 


inst 

N 


o 

ra 

t 

> 


m. circumferential location - radians 

* circumferential extent of distortion 
(largest region of pressure less than 
the average) 

* 3.14159 


* overall average 
« ring average 

=* applies to fan 

■ refers to time averaged value 

* applies to high pressure compressor 

* index of any ring counting frum hub 
to tip 

= refers to instantaneous average value 

* order of Fourier coefficient (number 
of vaves per circumference) 

■ reference 

* radial part applies to fan 

* total or stagnation 

* refers to engine inlet plane 
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